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SUMMARY 


-  "Several  years  of  usage  of  the  ARL  point  performance 
program  described  in  ARL  Technical  Engineering  Report  160  have 
indicated  the  desirability  of  certain  improvements  and  extensions 
which  have  been  incorporated.  This  document  describes  the  purpose 
and  nature  of  the  various  changes  in  sufficient  detail  to  allow 
it  to  be  used  independently  of  the  original  report. 
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Unit  or  Value 


Drag  coefficient 
Store  drag  coefficient 


Lift  coefficient 


Maximum  lift  coefficient 


Total  aerodynamic  drag 
Total  energy 


Energy  state 


Effective  thrust  vector 


Engine  net  thrust 
Engine  gross  thrust 
Acceleration  due  to  gravity 


Pressure  altitude 


Total  lift  force 


Jet  induced  lift  loss 


N(lb) 

J  (ft  lb) 
m  ( ft ) 

N(lb) 

N(  lb) 

N(lb) 

9.80665  m/s2 (32.17405  ft /s2  ) 
m(ft ) 

N(  lb) 

N(lb) 


Mach  number 

Load  factor  normal  to  aircraft 
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Specific  excess  power,  energy  m/s  (ft/s) 
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Intake  momentum  drag 
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N/m2  (lb /ft2  ) 
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Symbols 
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V 
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wf 


Definition 

Time 
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Fuel  flowrate 


Greek 
Symbol s 

a  Body  incidence  to  flight  path 

;  Load  factor  n  or  turn  rate  u 

according  to  selected  grid 
variable 

£'i  Engine  incidence  angle 


Nozzle  deflection  angle 


>,  Energy  rate  function 

t  Energy  state  Es  or  pressure 

altitude  hp  according  to 
selected  grid  variable 


Unit  or  Value 
s 

m/s(ft/s) 

N(  lb) 

kg/s( Ib/hr ) 

rad(deg) 

rad(deg) 

rad(deg) 

rad(deg) 
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Turn  rate 


deg/ s 
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1 .  INTRODUCTION 


Reference  1  described  a  suite  of  FORTRAN  programs  developed  for  the 
computation  and  presentation  of  data  used  in  evaluating  combat  aircraft 
performance. 

The  suite  of  programs  produces  tabulated  and  plotted  data  describing 
aircraft  combat  performance  in  energy  manoeuvrability  terms.  For  each  aircraft 
to  be  evaluated  the  user  must  supply  subroutines  which  evaluate  propulsion  and 
aerodynamic  characteristics. 

After  using  the  programs  for  several  years  some  improvements  and 
extensions  have  appeared  necessary  and  have  now  been  incorporated  in  the  suite 
of  programs. 

This  present  document  describes  the  purpose  and  nature  of  the  various 
changes  in  sufficient  detail  to  allow  it  to  be  used  independently  of 
Reference  1 . 

Chapter  2  presents  an  overview  of  the  suite  of  programs  and  a  description 
of  the  available  options  and  capabilities. 

Chapter  3  describes  the  modifications  to  the  mathematical  equations. 

Chapters  4,  5  and  6  present  a  self  contained  and  updated  user's  guide 
with  instructions  on  how  to  operate  the  three  modified  programs. 

Chapters  7,  8  and  9  describe  the  three  modified  programs;  details  of  the 
modifications  and  the  part  of  the  programs  they  affect  are  contained  in 
Appendix  A. 


Chapter  10  shows  examples  of  the  use  of  the  various  options. 


!.  OVERVIEW  AND  CAPABILITIES  OF  THE  SUITE  OF  PROGRAMS 

The  object  of  the  suite  of  programs  is  to  produce  tabular  and  plotted 
data  for  accurate  assessment  of  the  performance  and  manoeuvrability  of  any 
aircraft  operating  in  a  wide  range  of  weapon/store  configurations. 

The  suite  consists  of  three  main  programs  whose  interrelations  are  shown 
in  Figure  2.1. 

ACRAFT  is  the  generic  name  given  to  the  family  of  programs  Uiich 
represents  a  particular  aircraft;  it  includes  user  defined  routines  which 
calculate  propulsion  and  aerodynamic  parameters  using  engine  and  drag  data 
files.  The  program  uses  input  data  describing  the  aircraft  configuration  to 
perform  basic  calculations  and  produces  tabular  output  for  line  printer  listing 
as  well  as  a  data  file  for  input  to  program  P2V2. 

P2V2 ,  the  second  program  in  the  suite,  processes  the  data  provided  by 
program  ACRAFT.  For  some  requirements,  (options  4B,  4C ,  4Z)  plotter  files 
(P2.PLT)  are  produced  for  plotting  directly  on  the  off  line  plotter;  for  other 
requirements,  (4A,  4D,  4G)  program  P2V2  selects  and  rearranges  the  data  arrays 
for  further  processing  by  P4V2  . 

Program  P4V2  is  a  contour  plotting  program  which  processes  the  data 
provided  by  program  P2V2  and  produces  plotter  files  (P4.PLT)  for  submission  to 
the  off  line  plotter  queue. 

An  overview  of  all  the  options  available  is  shown  in  Figure  2.2. 

2 . 1  Program  ACRAFT 

Program  ACRAFT  produces  two  types  of  output: 

(1)  Unoptimised  grid  calculations: 

>  and  z,'  tabulated  over  an  unoptimised  grid  of  M,  4) 

(2)  Optimised  grid  calculations: 

X  t  and  C'  t  tabulated  over  an  optimised  grid  of  (Es.  4) 

where  the  symbols  can  represent  the  following  alternative  parameters 
is  energy  state  Es  or  pressure  altitude  hp. 
is  Mach  number. 

is  turn  rate  -  or  load  factor  n. 

is  load  factor  if  4  is  turn  rate, 
or  turn  rate  if  Z  is  load  factor. 

is  Ps,  Ps/wf ,  PsV/wf  or  Ps/V. 

is  the  maximum  value  of  >  over  the  permissable  range  of  M. 
is  the  z'  corresponding  to  XQpt 


M 

r 

C 

K 

X 

opt 

r  * 

'  Opt 
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The  selection  of  the  alternative  options  for  ■%,  t  and  £  is  made  through 
the  conversational  dialogue  when  running  the  program. 

In  addition  ACRAFT  provides  for  each  Mach  number  M  and  energy  state  x  the 
values  of  the  load  factor  and  turn  rate  corresponding  to  the  lift  limit 
boundary . 

2  .2  Program  P2V2 

The  input  to  the  program  P2V2  consist  of  data  files  prepared  for 
particular  aircraft  by  program  ACRAFT.  The  user  selects  by  terminal  input  the 
operations  to  be  performed  on  the  data.  The  user  can  choose  the  following 
options  for  presenting  the  data: 

(1)  Option  4.4:  Contours  of  x  on  a  (M,t)  grid  for  selected  ;  values. 

(2)  Option  4B:  Plots  of  X  versus  b  for  lines  of  constant  M  with  selected  t 

values. 

(3)  Option  4C:  Plots  of  Xopt  versus  Es  for  lines  of  constant  ; .  For 
this  option  the  second  type  of  output  produced  by  ACRAFT  is  necessa 
(see  paragraph  2.1). 

(4)  Option  4D:  Differential  contours  of  \  (between  two  aircraft)  on  a 
(M,t)  grid  for  selected  r  values. 

(5)  Option  4E:  Differential  4C  plots.  This  option  is  not  currently 
available. 

(6)  Option  4G:  Contours  of  •,  on  a  (M,;)  grid  for  selected  values  of  T. 

(7)  Option  4Z:  Plots  of  x  versus  M  and  versus  M  for  lines  of  constant  C 

with  selected  t  values. 

Note:  Now  that  option  4G  has  been  added  to  the  program,  option  4Z, 

(which  was  previously  used  to  create  manually  the  contours  of  4G),  is 
not  very  meaningful. 

When  selecting  option  4B,  4C  or  4Z,  the  output  from  the  program  is  a 
plotter  file  P2.PLT  which  can  be  submitted  to  the  off  line  plotter.  The  outputs 
produced  when  selecting  the  other  options  need  to  be  further  processed  by  the 
contour  plotting  program  P4V2. 

2 . 3  Program  P4V2 

Program  P4V2  processes  files  created  by  program  P2V2  (with  the 
extension  .CON)  when  selecting  options  4A ,  4D  and  4G  and  produces  a  plotter  file 
P4.PLT  to  be  submitted  for  off  line  plotting.  The  contour  plots  on  file  P4.PLT 
are  of  two  types: 

(1)  Contours  of  energy  rate  (or  other  energy  parameter)  on  a  Mach  number, 
altitude  (or  energy  state)  grid.  (Option  4A  and  4D) 

(2)  Contours  of  energy  rate  (or  other  energy  parameter)  on  a  Mach  number, 
turn  rate  (or  load  factor)  grid.  (Option  4G) 

A  conversational  dialogue  allows  the  user  to  define  the  range  of  contours 
required  and  the  scales  to  be  used  on  the  x  and  y  axes. 


[3) 


f 


OPTION 


4A  contour  of  \  on  a  (M,hp)  grid.  One  page  per  selected  n 


4B  x  versus  i>  for  lines  of  constant  M.  One  page  per  selected  2 


4D  differential  4A 


4G  contours  of  x  on  a  (M,n)  grid.  One  page  per  selected 


for  lines  of  constant  n.  One  page  per  selected  hp 


4A  ’  contours  of  \  on  a  (M,hp)  grid.  One  page  per  selected  * 


4B  x  versus  n  for  lines  of  constant  M.  One  page  per  selected  7 


1 

4D  differential  of  4A 


:  X  on  a  (M,i)  grid.  One  page  per  selected  9 


for  lines  of  constant  u. .  One  p 


4C  ;  x  opt  versus  E_  for  lines  of  constant  n 


differential  of  4C  -  Not  available 


4C  j  v  opt  versus  Eg  for  line  of  constant 


4E  differential  of  4C  -  Not  available 


contours  of  \  on  a  (M,E  )  grid.  One  page  per  selected  n 


4B  |  x  versus  *  for  lines  of  constant  M.  One  page  per  selectea  E 
1  s 

i 


4D  differential  4A 


G  contours  of  x  on  a  (M,n)  grid.  One  page  per  selected  E 


X  versus  M~)  for  lines  of  constant  n.  One  page  per  selected  Es!19 
*  versus  M  l  \ 


4A  I  contours  of  \  on  a  (M,Eg)  grid.  One  page  per  selected  1 


4B  x  versus  n  for  lines  of  constant  M.  One  page  per  selected  Es 


4D  differential  of  4A 


4G  I  contours  of  x  on  a  (M,„)  grid.  One  page  per  selected  Es 


X  versus  M 
n  versus  M 


for  lines  of  constant  One  page  per  selected  E 


'  is  energy  state  Es  or  pressure  altitude  hp 
£  is  load  factor  n  or  turn  rate  1 

K  is  energy  parameter  Pg  or  Pg/wf  or  Ps  V/wf  or  Ps/V 

FIG.  2.2:  Overview  of  Options  Available 


3.  ALTERATIONS  TO  EQUATIONS 


The  equations  of  motion  contained  in  Chapter  2  of  Reference  1  were 
approximate  in  that  net  thrust  FN  was  used  instead  of  gross  thrust  Fg  and  a 
ram  drag  term  was  not  included.  In  addition  the  effects  of  vectoring  the  thrust 
(as  occurs  in  STOVL  aircraft)  had  not  been  accounted  for;  this  involves  the 
inclusion  of  the  nozzle  deflection  angle  ©j  and  a  term  £L  to  allow  for  jet 
induced  lift  loss . 

The  equations  have  now  been  modified;  the  aerodynamic  and  propulsive 
forces  acting  on  the  aircraft  have  been  resolved  into  three  vectors  (see 
Figure  3.1): 

.  A  total  lift  vector  which  is  perpendicular  to  the  flight  path: 

Ltot  =  qSCL  +  (Fn  +  Rd)  sin  (a  +  0^  •+  i  j )  -  JL 


.  An  effective  thrust  vector  which  is  the  component  of  thrust  acting  along 
the  flight  path  (including  a  propulsive  drag  term): 

fefft  =  <fN  +  RD>  cos  +  ~'i  +  ©j)  ~  rd 

.  A  drag  vector  which  lies  along  the  flight  path: 

°TOT  =  qSCD  +  q  -  'jCDS 

=  Engine  net  propulsive  force 

*  Ram  drag  (Intake  momentum  drag) 

=  FN  +  RD  =  Engine  total  gross  thrust 

*  Angle  of  attack  (deg) 

*  Engine  incidence  angle  (deg) 

=  Nozzle  deflection  angle  (deg) 

”  Dynamic  pressure 

*  Jet  induced  lift  loss 

*  Lift  coefficient 

*  Drag  coefficient 


) 


where  F 


N 

Rd 

fg 

a 

3i 

?j 

q 

JL 

CL 


f 


We  can  then  write  as  with  conventional  aircraft  (see  Reference 


V 

(FEFFT  -  dtot)  - 
w 


(7] 


4.  PROGRAM  ACRAFT  USER'S  GUIDE 

The  conversational  input  to  the  program  is  best  described  with  reference 
to  a  typical  run.  Examples  of  dialogue  are  shown  in  Chapter  10. 

Care  must  be  taken  for  the  replies  to  questions  on  selecting  grid 
variables  to  be  compatible  with  the  format  of  the  plots  required  from  the 
subsequent  programs;  the  various  options  available  are  listed  in  Figure  2.2. 

Although  only  some  of  the  questions  have  been  modified  from  the  ones 
described  in  Reference  1,  for  completeness  all  replies  have  been  described  in 
this  updated  user's  guide. 

Question  1  :  IMPERIAL  (0)  or  SI(1)  UNITS 

Replies  of  0  or  1  indicate  that  output  quantities  from  ACRAFT  are  to  be 
in  Imperial  or  SI  units  respectively.  Any  other  reply  results  in  an  error 
message  and  the  question  is  repeated. 

Note:  All  internal  calculations  are  in  imperial  units. 

Question  2  :  AIRCRAFT  DATA  FILENAME 

Reply  with  the  name  (up  to  10  characters)  of  the  disk  file  containing  the 
required  data.  A  typical  example  of  a  store  file  is  given  in  Figure  4.1.  The 
contents  and  format  of  the  file  are  as  follows. 


Line  Format 


Descript  ion 


1  2A5 

2  3G 


3  3G 


4  14A5 

5  2A5 

6  2A5 

7  G 

8  10G 

9  10G 


Aircraft  name  (up  to  10  characters) 

Wing  reference  area  (ft  or  m  ) , 

Store  file  unit  flag  (0  for  Imperial  units,  1  for 
SI  units) 

Data  type  flag  (0  for  binary  spline  data,  1  for  PACAM 
tabular  data) 

Gross  weight  (lb  or  kg) 
c.g.  position  (2  MAC) 

Structural  load  factor  limit  of  aircraft 
Configuration  description  (up  to  70  characters) 

Engine  data  file  name  if  spline  data  or  PACAM  input 
file  name  if  tabular  data  (up  to  10  characters) 
Aerodynamic  data  file  name  if  spline  data  and  blank  if 
PACAM  tabular  data  (up  to  10  characters). 

Number  of  points  in  store  drag  table 
Mach  Number  list  in  store  drag  table  * 

Drag  count  list  in  store  drag  table  * 


*  10  items  per  line;  carry  on  to  next  line  if  more  than  10  entries  on  any 
line.  Not  needed  if  0  points  indicated  at  line  7. 
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Question  3  :  PS,  PS/WF,  PS*V/WF  OR  PS/V  (1,2,3  or  4) 


Reply  Meaning 

1  Calculate  Ps  as  dependent  variable,  i.e.  X  =  Ps 

2  Calculate  Ps/vf  as  dependent  variable,  i.e.  x  =  Ps/wf 

3  Calculate  PgV/1000wf  as  dependent  variable,  i.e.  x  =  psv/wf 

4  Calculate  1000PS/V  as  dependent  variable,  i.e.  x  =  Pg/V 

Values  of  1,2,3  and  4  are  valid  as  input.  Any  other  value  produces  an 

error  message  and  the  question  is  repeated. 

Question  4  :  IS  HP  THE  HEIGHT  VARIABLE 

A  reply  "Y"  (yes)  indicates  that  pressure  altitude  is  to  be  used  as  the 
height  grid  variable  (i.e.  t  is  pressure  altitude  hp).  Any  other  reply  will 
result  in  energy  state  being  used  as  the  height  variable  (i.e.  T  is  Es), 

If  energy  state  is  selected  as  the  height  variable  ( t  is  Es)  the 
following  question  is  asked. 

Question  4a  :  MAXIMUM  MANOEUVRE? 

A  reply  "Y"  will  cause  X  optimum  to  be  calculated  for  a  maximum  manoeuvre 
diagram  (MMD) . 

Any  other  reply  causes  the  performance  variable  x  to  be  calculated  for  an 
unoptimised  grid. 

Question  5  :  PRESET  GRID  (Y,N  or  C/R) 

A  Carriage  return  (C/R)  reply  results  in  the  preset  grid  values  being 
typed  on  the  user's  terminal.  The  values  for  the  preset  grid  are: 

start  value  increment  final  value 

t  (Imperial)  0  4000  68000 

’  (SI)  0  1000  20000 

M  (Military  Power)  0  .05  1.5 

M  (After  Burner)  0  0.05  2.0 

t  is  load  factor  1  2  9 

A  reply  "Y"  results  in  the  preset  grid  values  to  be  used  in  the  run.  The 
next  question  then  being  question  10. 

Any  other  reply  will  result  in  the  following  questions 
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Question  6  :  HP  ,HP  STEP,  NO  OF  STEPS  (if  T  is  h^) 

or  ES  ,ES  STEP,  NO  OF  STEPS  (if  T  is  Eg) 

Reply  with  .  an  initial  value  for  i  (must  be  >  0) 

.  an  increment  for  t  (must  be  >0) 

.  number  of  points,  including  first  and  last,  (must  be  >0  and 
<  50) 

These  three  parameters,  defining  the  height  grid,  should  be  typed  on  one 
line  separated  by  commas  or  blanks. 

Question  7  :  MACH,  MACH  STEP,  NO  OF  STEPS 


Reply  in  a  similar  way  to  question  6  so  as  to  specify  the  Mach  number 


If  an  optimized  grid  (MMD)  has  been  requested  in  question  4a,  question  7 
is  not  asked. 

Question  8  :  GRID  FOR  LOAD  FACTOR  (0)  OR  TURN  RATE  (1) 

A  reply  of  0,  requests  a  constant  increment  grid  of  load  factor  (C  is 
load  factor) 

A  reply  of  1,  indicates  a  request  for  a  constant  increment  grid  of  turn 
rate  (c,  is  turn  rate). 

Question  9  :  TR,  TR  STEP,  NO  OF  STEPS  (if  ;  is  turn  rate) 

or  GN0,  GN  STEP,  NO  OF  STEPS  (if  ?  is  load  factor) 

Reply  in  a  similar  fashion  to  question  6  so  as  to  specify  the  l  (turn 
rate  or  load  factor)  grid. 

Note:  Tumrate  is  expressed  in  deg/s  and  if  5  is  load  factor  the  initial  value 

must  be  greater  or  equal  to  1. 

Question  10  :  POWER  (MIL-100,  MAX-200) 

Military  power  is  indicated  by  a  value  of  100,  and  maximum  augmented 
power  by  a  value  of  200.  Tens  and  units  digits  indicate  a  percent  of  the 

maximum  of  that  range.  A  valid  reply  is  0  <  Power  <  200.  Any  other  reply 

causes  an  error  message  to  be  typed  on  the  terminal  and  the  question  is 
repeated. 

Note:  The  power  setting  requested  must  be  consistent  with  the  propulsion 

character ist ics  specified  in  the  user  defined  subroutine  THRUST. 

Question  11  :  WING  SWEEP  (IF  VARIABLE) 


Supply  wing  sweep  in  degrees  if  aerodynamic  data  includes  wing  sweep 
effects.  Otherwise  reply  with  a  carriage  return. 

A  nominal  valid  reply  is  0  <  Wing  sweep  <  80*.  Any  other  reply  causes  an 
error  message  to  be  typed  on  the  terminal  and  the  question  is  repeated. 
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Due st ion  12  :  NOZZLE  DEFLECTION  (IF  VARIABLE) 


Supply  the  nozzle  deflection  in  degrees  if  the  aircraft  is  a  vectored 
thrust  aircraft;  otherwise  reply  with  a  carriage  return. 

A  nominal  valid  reply  is  0  <  nozzle  angle  <  110*.  Any  other  reply  causes 
an  error  message  to  be  typed  on  the  terminal  and  the  question  is  repeated. 

Question  13  :  ATMOSPHERE,  DEVIATION 

This  input  is  provided  if  the  user-defined  subroutine  THRUST  (from 
ACRAFT)  has  been  written  to  handle  alternative  atmospheres.  Valid  replies  are: 

ICAO  ,X  *1  format  (A4,G) 

ARDU.X  J 

Where  X  is  a  deviation  from  the  nominated  atmosphere  in  degrees  Celsius. 
It  may  be  omitted  if  X  is  zero. 

If  an  invalid  reply  is  made  a  warning  message  is  typed  to  the  terminal 
and  ICAO  atmosphere  is  assumed. 


Question  14  :  OUTPUT  (1*TEXT,  2=NOS,  3-B0TH) 


Reply  Meaning 

1  Formatted  output  with  full  ASCII  text  and  headings 
is  requested  (on  logical  unit  6). 

2  Numerical  data  output  is  required  (on  logical  unit 
8)  for  input  to  program  P2V2. 

3  Both  types  of  output  are  required. 

Any  other  reply  causes  an  error  message  to  be  typed  and  the  question  is 
repeated . 

Question  15  :  0/P  FILE  NAME  FOR  UNIT  m 


Supply  a  10-character  filename  for  output  as  requested;  m  *  6  is  the 
logical  unit  number  for  full  format  output;  m  *  8  is  the  logical  unit  number  for 
numerical  output . 

Formatted  output  should  be  printed  with  /P/B  switches  to  produce  listings 
suitable  for  permanent  retention. 

Numerical  output  may  be  printed  if  required. 


! 
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AIRCRAFT 

230. . 0.0 

25000. . 00. .7.0 

AIR  SUPERIORITY  :  2  AIM9  +  50?  FUEL 
ENGIffE.BIN 
AERO.BIN 
5 

•6, .7, .8, .9, 1.0 

120,1*6,191, 304,370 


FIG.  4.1  Typical  Example  of  a  Store  File 


( 
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5.  PROGRAM  P2V2  USER’S  GUIDE 


The  inputs  to  program  P2V2  consist  of  the  data  files  created  by  the 
program  ACRAFT. 

The  data  is  then  manipulated  by  P2V2  according  to  commands  supplied  at 
the  users  terminal  during  the  execution  of  the  program. 

Question  1  :  OPTION  OR  (CR)  FOR  HELP: 

Reply  with  the  appropriate  option  code.  A  carriage  return  is  interpreted 
as  a  request  for  help  and  a  one  line  description  of  each  option  together  with 
the  name  of  the  output  file  produced,  is  typed  on  the  terminal. 

Reply  Option 

4A  X  contours  on  a  (M,  t)  grid  for  selected  ;  values. 

4B  X  versus  4'  for  lines  of  constant  M  for  selected  t  values. 

4C  Xopt  versus  Es  f°r  constant  4  values. 

4D  differential  x  contours  on  a  (M,t)  grid  for  selected  4  values. 

4E  differential  \Q  versus  Es  for  constant  4  values.  Not 

implemented  in  this  version. 

4F  Maximum  Manoeuvre  Persistence  interactive  aid.  Not  available 

in  this  version. 

4G  X  contours  on  a  (M, 4)  grid  for  selected  '  values. 

4Z  X  versus  M  for  lines  of  constant  4  as  well  as  4'  versus  M  for 

lines  of  constant  4  for  selected  t  values. 

any  other  reply  will  cause  a  help  message  to  be  typed  to  the 
terminal  and  question  1  to  be  repeated. 


Where  X  is  Pg  or  Pg/wf  or  PgV/wf  or  Ps/V. 
T  is  Eg  or  pressure  height  hp. 

4  is  turn  rate  or  load  factor. 


4'  is  load  factor  if  4  is  turn  rate  or  vice  versa. 
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5.  PROGRAM  P2V2  USER'S  GUIDE 


The  inputs  to  program  P2V2  consist  of  the  data  files  created  by  the 
program  ACRAFT. 

The  data  is  then  manipulated  by  P2V2  according  to  commands  supplied  at 
the  users  terminal  during  the  execution  of  the  program. 

Question  1  :  OPTION  OR  (CR)  FOR  HELP: 

Reply  with  the  appropriate  option  code.  A  carriage  return  is  interpreted 
as  a  request  for  help  and  a  one  line  description  of  each  option  together  with 
the  name  of  the  output  file  produced,  is  typed  on  the  terminal. 


Reply 

Opt 

ion 

4A 

\  contours  on  a 

(M, 

t  )  gr i d 

for  selected  : 

;  values  . 

4B 

X  versus  £ '  for 

lines  of  con 

st  ant 

M  for  se 

lected  i  values 

• 

4C 

vopt  versus  Es 

for 

constant 

l  va 

lues. 

4D 

differential  v 

contours  on  a 

(M,t 

)  grid  for  selected  l  va 

1  ues 

4E 

differential  y 

versus  Es 

for 

constant 

C  val ues  .  Not 

implemented  in 

fh  is 

vers i on . 

4F 

Maximum  Manoeuvre  Persistenc 

e  int 

erac  t i ve 

aid.  Not  avail 

able 

in  this  version 

• 

4G 

.  contours  on  a 

(M.C 

)  grid  for  se 

lected  ~ 

val ues  . 

4Z 

*  versus  M  for 

1  i  nes 

of  cons 

t  ant 

s  as  well 

as  versus  M 

for 

lines  of  constant  r 

for  selected 

t  values. 

- 

any  other  reply 

wi  1 1 

cause  a 

help 

message 

to  be  typed  to 

the 

terminal  and  question  1  to  be  repeated. 
Where  <  is  Pg  or  Pg/wf  or  PgV/wf  or  P$/V. 

*  is  Es  or  pressure  height  hp. 

^  is  turn  rate  or  load  factor. 

~ '  is  load  factor  if  ^  is  turn  rate  or  vice  versa. 
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5 . 1  Option  4A 

This  option  produces  X  contour  data  on  a  Mach  number  (x  axis),  t  (y  axis) 
grid,  at  selected  values  of  ;  (turn  rate  or  load  factor).  The  output  is  written 
on  file  P24A1 .CON  which  will  have  to  be  processed  by  P4V2  . 

Examples  of  dialogues  running  this  option  are  shown  in  Figures 

10.1  and  10.4. 

Question  Al  :  DATA  BASE  FILENAME 

Reply  with  the  name  (up  to  10  characters)  of  the  appropriate  numerical 
output  (logical  unit  8)  file  from  P1V2.  After  reading  the  header  on  the  file 
program  P2V2  echoes  on  the  user  terminal. 

***  DATA  ARE  IN  V  I'NITS,  ENERGY  PARAMETER  IS  \  *** 

***  CONSTANT  INC.  OF  l  *** 

U  is  IMPERIAL  or  SI ,  \  is  PS,  PS/WF,  PS.V/WF  or  PS /V  as  appropriate  and  Z 
is  Tl’RNRATE  or  LOAD  FACTOR  as  appropriate. 

Question  A2  :  REPLY  "YES",  "NO",  "ALL"  OR  "END" 

OUTPUT  DATA  FOR  ;  =  ^  ? 

Reply  of  ' Y ‘ ,  'N',  'A'  or  ’ E '  indicates  as  follows: 

Repl y  Meaning 

Y  Output  data  for  contour  plot  at  Z  ~Zi  required. 

N  No  output  data  required  for  Z  *  ^ 

A  Output  data  required  for  all  . 

E  Exit  from  output  data  loop  and  do  not  output  data  for 

C  =  Z  ^  or  any  following 

Any  other  reply  is  interpreted  as  a  'N'. 

Question  A2  is  repeated  for  all  values  of  t.  selected  unless  the  loop  is 

ex  ited . 


Output  is  written  onto  file  P24A1.CON,  for  input  to  program  P4V2 .  File 
P24A1 .CON  may  be  renamed  if  desired  before  running  P4V2.  P24A1 .CON  is  an  ASCII 

file  and  may  be  printed,  although  it  is  not  easy  to  read  but  is  suitable  for 
debugging  purposes. 
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5 . 2  Option  4B 

This  option  creates  plots  of  \  versus  C  '  for  selected  values  of  Mach 
number.  One  page  of  plots  is  produced  for  each  requested  altitude  (or  energy 
state).  The  output  is  written  onto  the  plot  file  P2.PLT  v*iich  can  be  submitted 
to  the  plotter  queue. 

An  example  of  dialogue  running  this  option  is  shown  in  Figure  10.7. 
Question  Bl  :  DATA  BASE  FILENAME: 

Reply  with  the  name  (up  to  10  characters)  of  t..e  appropriate  numerical 
output  (logical  unit  =  8)  file  from  P1V2.  Program  P2V2  will  echo  on  the 
terminal : 

***  DATA  ARE  IN  U  UNITS ,  ENERGY  PARAMETER  IS  '  *** 

***  CONSTANT  INCREMENT  OF  •'  *** 

where  (as  appropriate) 

C  is  IMPERIAL  or  SI 

is  PS,  PS/WF,  P  S . V  '  WF  or  PS/V 
is  TURN  RATE  or  LOAD  FACTOR 

Question  B2  :  SCALES  IN  UNITS/ INCH  OF  PLOT 

TURN RATE /LOAD  FACTOR,  ENERGY  VARIABLE: 

Reply  with  the  amount  of  physical  quantity  which  one  inch  of  plot  would 
represent  on  each  axis  where  the  x  axis  represents  Z'  (turnrate  or  load  factor) 
and  the  y  axis  represents  the  energy  variable  \  .  (As  a  rough  guide  use  500 
ft/s/inch  for  Ps  and  A  deg/sec/inch  for  the  turn  rate). 

Question  B3  :  MAXIMUM  FOR  ENERGY  AXIS: 

Reply  with  the  approximate  maximum  \  value  of  interest  on  the  y  axis  (in 
the  order  of  1000  ft/s  for  Ps). 

Question  BA  :  ENERGY  AXIS  LENGTH: 

Reply  with  the  appropriate  length  in  inches  of  the  y  axis  (7  inches  is 
representative  for  trimming  to  AA  size). 

Question  B5  :  REPLY  "YES",  "NO",  "ALL"  OR  "END": 

HEIGHT  =  Ti  ? 

One  page  of  curves  is  plotted  for  each  height  selected. 

Replies  and  significance  of  questions  are  listed  below. 
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Mean i ng 


Reply 

Y  (Yes) 

N  (No) 

A  (All) 

E  (End) 

Quest  ion  B6  :  MACH  =  Mj 

This  question  is  asked 
Reply 

Y  (Yes) 

N  (No) 

A  (All) 

E  (End) 

The  output  file  P2.PLT 
the  off  line  plot  queue. 


Plotting  is  required  for  value.  Mach  numbers 

to  be  plotted  will  now  be  listed. 

Do  not  produce  the  plot  for  this  value  of  t.. 

Plot  all  the  data  given  for  Use  this  when  the 

Mach  number  values  are  known  and  plots  for  all 
are  to  be  produced. 

Exit  from  height  plotting  and  hence  from  program. 

•> 

only  if  the  reply  to  Question  E5  was  ' Y'. 

Mean i ng 

Plot  a  curve  for  this  Mach  number. 

Skip  this  Mach  number. 

Plot  curves  for  all  remaining  Mach  numbers  for 
current  T;  . 

Skip  remaining  Mach  numbers  and  go  to  the  next  t,  . 
contains  the  reauested  plots  and  is  submitted  to 
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5 . 3  Option  4C 

This  option  prodaces  plots  of  v0pt  versus  energy  state  Es  for  content!, 
values.  The  output  is  written  onto  the  plot  file  P2.PLT  which  can  be  submitted 
to  the  plot  queue. 

Examples  of  dialogues  running  this  option  are  shown  in  Figures 
10.9  and  10.12. 

Question  Cl  :  DATA  BASE  FILENAME: 

Reply  with  the  name  (up  to  10  characters)  of  the  numerical  output  file 
from  P1V2.  Program  P2V2  will  echo  the  following  message  on  the  terminal: 

***  DATA  ARE  IN  U  UNITS,  ENERGY  PARAMETERS  IS  \  *** 

***  CONSTANT  INC.  OF  -  **** 

where  U  is  IMPERIAL  or  SI 

X  is  PS  or  PS/WF  or  PS.V/WF  or  PS/V 
’  is  TURN  RATE  or  LOAD  FACTOR 

Question  C2  :  SCALE  IN  UNIT  'INCH  OF  PLOT 
ES,  PS  : 

Reply  with  the  amount  of  physical  quantity  which  one  inch  of  plot  would 
represent  on  each  axis.  Es  is  plotted  along  the  x  axis  and  .  is  plotted  along 
the  y  axis. 

Question  C3  :  PS  MIN  VALUE,  PS  AXIS  LENGTH  (INS)  : 

Reply  with  the  approximate  minimum  »  value  of  interest  and  the  length  of 
the  v  axis  tin  inches). 

Question  C4  :  REPLY  "YES",  "NO",  "ALL"  OR  "END" 

z  =  * .  ? 

-  "1 

Where  !  is  TURN  RATE  or  LOAD  FACTOR  as  appropriate. 

Replies  and  meanings  are  indicated  below. 


Reply  Mean i ng 

Y  (Yes)  Plot  a  curve  of  (Es,  v)  points  for  this 

N  (No)  Do  not  plot  a  curve  and  try  next  . 

A  (All)  Plot  curves  for  this  and  subsequent  . 

E  (End)  Skip  the  remaining  and  hence  exit  from  program. 

The  output  file  for  option  4C  is  P2.PLT,  a  plot  file  which  must  be 

submitted  to  the  off-line  plot  queue. 
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5 .A  Option  AD 

This  option  produces  differentia!  i  contour  data  on  a  Mach  number,1 
grid  at  selected  values  of  :  (turn  rate  or  load  factor).  The  versions  of  ACRAFT 
appropriate  to  the  two  different  aircraft  are  run  using  their  respective  data 
files  with  identical  altitude,  Mach  number  and  load  factor  (or  turn  rate)  grids 
as  in  option  AA ,  and  then  subtracting  the  data  at  the  grid  points  when  running 
P2V2.  The  output  is  written  on  three  files  P2AA1.C0N,  P2AA2.CON  and  P2ADIF.C0N 
which  will  have  to  be  processed  by  PAV2  . 

An  example  of  dialogue  running  this  option  is  shown  in  Figure  10.15. 
Question  D1  :  DATA  BASE  FILENAME: 

Reply  with  the  name  (up  to  10  characters)  of  the  appropriate  numerical 
output  file  (logical  unit  S)  created  by  P1V2  for  the  first  aircraft.  After 
reading  the  header  on  the  file  program  P2V2  echoes,  as  in  option  AA,  details  of 
the  options  selected  when  running  P1V2. 

Question  D2  :  REPLY  "YES",  "NO",  "ALL"  OR  "END" 

OUTPUT  DATA  FOR  ;  =  ^7 

Replv  to  this  question  is  the  same  as  for  option  AA.  The  load  factor  (or 
turn  rate)  selection  applies  both  to  the  base  file  and  to  the  comparison  file 
considered.  This  question  is  repeated  for  all  values  of  selected  when  the 
loop  is  exited. 

Question  D3  :  COMPARISON  FILENAME: 

Reply  with  the  name  of  the  output  file  (logical  unit  8)  created  by  P1V2 
for  the  second  aircraft  (comparison  aircraft)  to  be  compared  with  the  first 
aircraft  (base  aircraft).  Data  files  are  considered  to  be  valid  for  comparison 
if  when  running  ACRAFT  the  grid  points  selected  are  identical,  and  the  type  of 
energy  parameter  and  units  chosen  are  the  same;  other  parameters  such  as  power 
setting,  wing  sweeps,  nozzle  deflection  and  atmosphere  type  are  not  checked. 

Output  for  this  option  is  in  ASCII  mode  on  three  files  P2AA1.C0N, 
P2AA2.CON,  P2ADIF.C0N  which  can  be  used  as  input  for  processing  by  PAV2 . 

P2AA1 .CON  is  the  same  as  would  be  produced  by  running  option  AA  with  the  same 
replies  using  the  base  aircraft  data  file  as  input.  P2AA2.CON  is  the  same  as 
would  be  produced  by  running  option  AA  using  the  comparison  aircraft  data  file 
as  input.  P2ADIF.C0N  is  the  file  used  when  running  PAV2  to  generate  the 
differential  contour  plots. 

5 . 5  Option  AE 


In  Reference  1  this  option  produced  differential  MMD  plots.  This  option 
was  not  implemented  when  the  programs  were  restructured  and  extended  as  it  was 
felt  the  amount  of  programing  effort  required  was  not  in  proportion  to  the 
usefulness  of  the  plots  presented.  The  code  used  in  Reference  1  is  still 
available  and  this  option  could  be  implemented  if  required  at  a  later  date. 
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5 .6  Option  4F 

In  Reference  i  this  option  provided  an  interactive  aid  for  maximum 
manoeuvre  persistence  where  the  user  provided  data  obtained  from  the  output  file 
of  option  4C  to  calculate  the  number  of  turns  an  aircraft  could  make  before 
running  out  of  fuel.  This  option  was  removed  because  it  was  felt  those 
calculations  were  better  handled  off  line  with  a  calculator. 

5.7  Option  4G 

This  option  creates  a  file  P24G.CON  which  can  then  be  processed  by 
program  P4V2  to  produce  contours  of  >,  on  a  Mach  number  (x  axis),  £  (y  axis)  grid 
at  selected  values  of  t. 

An  example  of  dialogue  running  this  option  is  shown  in  Figure  10.17. 
Question  Cl  :  DATA  BASE  FILENAME  : 

Reply  with  the  name  (up  to  10  characters)  of  the  appropriate  numerical 
output  file  from  P1V2.  After  reading  the  header  on  the  file,  program  P2V2  will 
echo  the  following: 

***  DATA  ARE  IN  U  UNITS,  ENERGY  PARAMETER  IS  *** 

***  CONSTANT  INC.  OF  ;  *** 

where  U  is  IMPERIAL  or  SI 

X  is  PS  or  PS/WF  or  PS.V/WF  or  PS/V 
£  is  TURN  RATE  or  LOAD  FACTOR 

Question  G2  :  REPLY  "YES",  "NO",  "ALL"  OR  "END" 

OUTPUT  DATA  FOR  T  -  T i  ? 


:  ' Y '  ,  ’N’  , 

'A',  or  'E*  indicates  as  follows 

Reply 

Meaning 

Y  (yes) 

Output  data  for  contour  plot  at  t  =  t^is 

requi red 

N  (no) 

No  output  data  required  for  t  = 

A  (all) 

Output  data  required  for  all 

E  (end) 

Exit  from  loop  and  program  P2V2. 

Do  not  output  data  for  or  any  subsequent 

•  i  r  .  i  i  _ i  f  _  j  _ 4.1 _ 

U  • 

1  _ _ 

Question  G2  is  repeated  for  all  values  of  t  selected  unless  the  loop 
is  exited. 

The  output  is  written  onto  the  file  P24G.CON  for  input  to  program  P4V2 . 
P24G.CON  is  an  ASCII  file  and  maybe  printed,  although  it  is  not  easy  to  read  but 
suitable  for  debugging  purposes. 
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5 .8  Option  42 

This  option  creates  a  plot  file  P2.PLT  for  plotting  the  energv  parameter 
X  versus  Mach  number  for  lines  of  constant  as  well  as  versus  Mach  number 
for  lines  of  constant  I  .  One  page  of  plots  is  produced  for  each  requested  t  . 

An  example  of  dialogue,  running  this  option  is  shown  in  Figure  10.21. 

Question  Z1  :  DATA  BASE  FILENAME: 

Reply  with  the  name  of  the  appropriate  numerical  output  file  from  P 1 V 2 . 
Program  P2V2  echoes  the  following  message  on  the  terminal: 

***  DATA  ARE  IN  U  WITS,  ENERGY  PARAMETER  IS  .  *** 

***  CONSTANT  INC.  OF  :  *** 

where  U  ,  \  and  '  are  appropriately  defined  as  for  option  4B. 

Question  Z2  :  REPLY  "YES",  "NO",  "ALL"  OR  "END" 

HEIGHT  *  ti  ? 

One  page  of  curves  is  plotted  for  each  height  selected.  Reply  as  for 
question  B5 

Question  23  :  C  =  ;  ? 

This  question  is  asked  only  if  the  answer  to  question  Z2  is  yes. 

Reply  as  for  question  B6  with  Mach  number  being  appropriately  replaced  bv 
in  the  discussion. 

The  output  file  for  option  4Z  is  P2.PLT,  a  plot  file  which  must  be 
submitted  to  the  off  line  plot  queue.  The  scales  for  this  option  are  fixed. 
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6.  PROGRAM  P4V2  USER'S  GUIDE 


The  inputs  to  program  P4V2  consist  of  data  files  created  by  program  P2V2 
when  choosing  options  4A ,  4D  and  4G. 

A  selection  is  made  when  running  P4V'2  using  terminal  inputs  of  data  to  be 
plotted.  Numerous  examples  are  shown  in  Chapter  10. 

The  sole  output  of  the  program  P4V2  is  the  plotter  file  P4.PLT. 

6 . 1  Option  4A 

Examples  of  dialogues  running  this  option  are  shown  in  Figures 

10.1  and  10.4. 

Question  1  :  INPUT  FILENAME 

Reply  with  the  name  (up  to  10  characters)  of  the  appropriate  file 

generated  by  option  4A  in  program  P2V2.  The  file  P24A1.C0N  created  by  P2V2  may 

be  renamed  before  running  P4V2 ,  (P4V2  will  still  know  which  option  produced  it) 
P4V2  will  then  echo  on  the  users  terminal  the  following  options  selected  in 
P1V2: 

***  DATA  ARE  IN  U  UNITS ,  ENERGY  PARAMETER  IS  \  *** 

***  CONSTANT  INC.  FOR  I  *** 

Where  (as  appropriate) 

U  is  IMPERIAL  or  SI 
X  is  PS,  PS/WF,  PS .V /WF  or  PS/V 
4  is  LOAD  FACTOR  or  TURN  RATE 

Question  2  :  SCALES  IN  UNITS/ IN  OF  PLOT  -  MACH,  ALT: 

Reply  with  the  amount  of  physical  quantity  which  one  inch  of  plot  should 

represent  on  each  axis,  where  the  x  axis  represents  Mach  number  and  the  y  axis 

represents  pressure  altitude  or  energy  state.  (As  a  rough  guide  use  .2  to  .4 
units/inch  for  the  Mach  number  and  10000  ft/inch  for  the  altitude). 

Question  3  :  SMOOTHED  AND  TEXTURED  CONTOURS: 

Reply  with  "Y"  for  the  contour  points  to  be  joined  by  a  cubic  spline, 
meaning  contours  will  have  continuous  1st  and  2nd  derivatives. 

Reply  with  "N"  for  the  contour  points  to  be  joined  by  straight  lines. 

If  the  grid  size  is  small  there  is  little  or  no  difference  between 
smoothing  and  non-smoothing  (So  it  is  better  to  select  non-smoothing  as  this 
will  speed  up  execution  time). 

Question  4  :  PLOT  ES  CONTOURS? 

Reply  with  "Y"  (yes)  or  "N"  (no)  where  yes  signifies  that  dashed  lines 
of  constant  energy  state  will  be  superimposed  on  the  requested  plots. 
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Question  5  :  CONTOUR  LEVELS  -  START,  STEP,  NO: 

This  question  is  only  asked  if  energy  state  contours  are  required. 

Reply  with  the  minimum  Eg,  the  increment  in  Es  between  contours  and 
the  number  of  Es  contours  required. 

These  values  should  be  in  the  appropriate  units  (as  a  rough  guide  use 
5000ft  for  the  start  and  step  values  and  20  for  the  number  of  steps). 

Question  6  :  PLOT  L  - 

Z  is  either  LOAD  FACTOR  or  TURN  RATE  as  appropriate. 

Reply  with  "Y"  (yes)  or  "N"  (no)  where  yes  signifies  that  a  new  page  of 
plots  be  set  up  for  the  required  J  .  If  the  reply  is  no,  the  question  is 
repeated  for  the  next 

Question  7  :  CONTOUR  LEVELS  -  START,  STEP,  NO: 

If  the  reply  to  question  6  was  "yes"  then  this  question  is  asked  for  the 
page  of  plot  for  each  requested  t,. 

Reply  with  the  \  value  of  the  lowest  valued  contour  to  be  plotted,  the 
increment  between  contours  (must  be  positive)  and  the  number  of  contours  (>0  and 
<  50). 


These  values  should  be  in  the  appropriate  units  depending  upon 
both  U  and  '<  .  As  load  factor  increases,  more  negative  contours  should  be 
plotted . 

After  the  last  value  of  s  is  processed,  an  informative  message  "END  OF 
DATA  ON  LOG5"  is  typed  on  the  terminal,  it  does  not  indicate  an  error 
cond  it  ion . 

After  question  7  is  repeated  for  each  Z  value  for  which  data  exists,  PiV2 
ceases  execution  and  the  output  plot  file  P4.PLT  may  be  placed  on  the  off  line 
plotter  queue. 

6 .2  Option  4P 


An  example  of  dialogue  running  this  option  is  shown  in  Figure  10.15. 

Replies  for  this  option  are  very  similar  to  those  of  option  4A  for  one 

aircraft.  The  only  difference  is  that  requests  for  contour  levels  now  refer  to 

the  difference  between  the  two  aircraft  at  each  grid  point.  The  aircraft 

referrred  to  by  program  P2V2  as  the  DATA  BASE  will  be  indicated  by  positive 
contours  and  that  referred  to  as  the  COMPARISON  will  be  indicated  by  negative 
contours.  Full  identification  headers  are  included  on  the  plotted  output  on 
file  P4.PLT. 


r 
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6 . 3  Option  4G 

An  example  of  dialogue  using  this  option  is  shown  in  Figure  10.17. 
Question  1  :  INPUT  FILENAME 

The  name  of  the  file  created  by  program  P2V2  when  using  option  4G  is 
P24G.C0N  but  this  file  can  be  renamed  before  running  P4V2 .  See  Paragraph  6.1 
for  the  meaning  of  the  prompt  message. 

Question  2  :  SCALES  IN  UNITS/ IN  OF  PLOT  -  MACH,  TR/LF: 

Reply  with  the  amount  of  physical  quantity  which  one  inch  of  plot  should 
represent  on  each  axis,  where  the  x  axis  represents  Mach  number  and  the  y  axis 
represents  either  turn  rate  or  load  factor.  (As  a  rough  guide  use  .2  to  .4 
unit/inch  for  the  Mach  number,  and  4  deg/sec/inch  for  the  turn  rate) 

Question  3  :  SMOOTHED  AND  TEXTURED  CONTOURS: 

Reply  with  ' Y 1  or  'N'  as  for  question  3  of  option  4A. 

Question  4  :  PLOT  LOAD  FAC.  /TURN  RAD.  GRID  ? 

This  question  is  only  asked  if  Z  is  turnrate. 

Reply  with  ' Y 1  (yes)  if  a  background  grid  of  constant  load  factor  lines 
and  constant  turn  radius  lines  is  desired  for  each  page  of  plot.  The  grid  is 
appropriate! v  modified  to  suit  the  unit  system  selected. 

Any  other  reply  will  be  interpreted  as  no. 

Question  5  :  PLOT  H  =  Ti? 

Reply  with  "Y"  (yes)  if  a  new  page  of  contours  is  to  be  set  up  for  the 
required  height. 

Any  other  reply  is  interpreted  as  no  and  the  question  is  repeated  for  the 
next  height  value. 

Question  6  :  CONTOURS  LEVELS  -  START,  STEP,  NO: 

If  the  reply  to  question  5  was  "yes"  then  this  question  is  asked  for  the 
page  of  plot  for  each  required  height. 

Reply  as  for  question  7  of  option  4A. 

After  question  6  is  repeated  for  each  height  value  for  which  data  exists, 
P4V2  ceases  execution  and  the  output  plot  file  P4.PLT  can  be  placed  on  the  off 
line  plot  queue. 
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7.  PROGRAM  ACRAFT  DESCRIPTION 


7 . 1  Loading  of  Program  on  DEC  System-10 

The  commands  for  loading  program  ACRAFT  and  saving  a  core  image  on  disk 
on  DEC  svstem-10  are: 

LOAD  ACRAFT,  P1V2,  PI  LIB./ SEARCH  ,  PAC2  IN/SEARCH 
SAVE  ACRAFT 

The  SAVE  command  on  the  DEC  system- 10  creates  an  executable  version  of 
the  program. 

The  file  ACRAFT  contains  all  the  user  defined  routines  which  are 
different  for  each  aircraft,  whereas  the  service  routines  in  P1V2  are  general 
and  applicable  to  all  aircraft.  P1LIB  includes  routines  for  use  with  spline 
fitted  data  (Reference  3)  whereas  PAC2IN  contains  routines  for  use  with  the 
tabular  format  data  of  a  PACAM2  aircraft  data  file  (Reference  2). 

The  routines  contained  in  each  of  these  files  are  listed  in  the  following 
table  and  described  in  this  Chapter. 


ACRAFT 

MAIN  BLOCK  THRUST  AERO MX  GETALF  AERO  GETCD 
GETCDCL  XSPMIN  XSPMAX  YSPMIN  YSPMAX 

GETCDA 

P1V2 

PI  IN 

PI  OUT 

BADINP 

BININ 

TABLE 

ALTIT 

ATMOS 

P ARAMS 

INTRP 

MAX  MAN 

HEIGHT  HTRUE 

SEP  P10UTA 

IDENT 

MON SEP 

PI  LIB 

SURF 

CUBICS 

CHECKD 

SPDER3 

ROM  IN 

PAC2IN 

PAC2IN 

TABIN 

PACINT 

FIND 

DINTRP 

PACDER 

7.2  File  ACRAFT 

7.2.1  ACRAFT  MAIN  Program 


This  is  a  small  routine  which  calls  P1IN  for  program  input  and  TABLE  tc 
start  calculation  of  aircraft  performance  parameters. 

7.2.2  BLOCK  DATA 

Defines  the  data  in  the  COMMON  blocks.  The  data  in  COMMON  C  represents 
atmospheric  profiles,  text  constants  and  other  numeric  constants;  the  data  in 
COJflON  USER  is  defined  by  the  user  for  use  in  the  aircraft  dependent  routines. 
The  contents  of  all  the  COMMON  blocks  are  listed  in  Appendix  6  of  Reference  1. 
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7.2.3  THRUST  Subroutine 


THRUST  is  an  aircraft  dependent  routine  supplied  by  the  user  and  is 
called  by  soubroutine  TABLE  or  SEP  (in  P1V2). 

In  THRUST,  values  of  gross  thrust  (and  lift  loss  in  the  case  of  vectored 
thrust  aircraft)  are  determined  and  transmitted  through  COMMON  USER  to  the  user 
defined  subroutines  AERO  and  AEROMX. 

Fuel  flow  and  maximum  Mach  number  are  also  calculated  in  THRUST  and 
supplied  to  P1V2  through  COMMON  B.  Since  THRUST  is  user  defined  its  form  will 
depend  on  the  nature  of  the  data  available.  The  following  three  examples 
demonstrate  the  flexibility  of  this  approach. 

In  example  1  (Figure  7.1)  spline  coefficients  have  been  derived  from 
original  data  for  net  thrust,  ram  drag  and  fuel  flow  expressed  as  dimensional 
functions  of  height  and  Mach  number.  The  data  was  available  for  two  power 
settings  (intermediate  and  augmented)  and  an  ICAO  atmosphere.  The  spline 
coefficients  were  derived  using  program  SPLFIT.EXE  (see  Reference  3)  and  stored 
in  the  binary  file  ENGIN.BIN  which  is  then  accessed  in  THRUST  by  calls  to  SURF. 

In  example  2  (Figure  7.2)  spline  fitting  techniques  were  applied  to 
tables  of  gross  thrust,  net  thrust  and  fuel  flow  data,  available  in  normalised 
form  as  functions  of  height  and  Mach  number;  an  additional  parameter,  lift  loss 
coefficient,  was  expressed  as  a  function  of  jet  velocity  ratio  and  nozzle 
deflection  angle. 

Example  3  (Figure  7.3)  illustrates  the  case  where  the  data  is  accessed 
from  an  ASCII  file  with  tabular  format.  Using  the  routines  from  PAC2IN,  gross 
thrust,  ram  drag  and  fuel  flow  are  linearly  interpolated  for  the  desired 
altitude,  Mach  number  and  power  setting. 

7.2.4  AEROMX  Subroutine 


AEROMX  is  a  user  defined  routine  called  by  subroutine  TABLE  or  SEP. 

For  a  given  Mach  number,  AEROMX  returns,  to  P1V2,  through  COMMON  B,  the 
maximum  usable  lift  coefficient  together  with  the  corresponding  angle  of  attack 
and  total  lift. 

Other  aerodynamic  information  dependent  only  on  height  and  Mach  number 
can  also  be  determined  in  AEROMX  (e.g.  incremental  store  drag  and  other 
aerodynamic  coefficients)  and  transmitted  to  the  other  user  defined  routines 
through  COWON  USER.  An  example  is  shown  in  Figure  7.4. 

7.2.5  GETALF  Subroutine 


GETALF  is  a  user  defined  routine  called  by  AEROMX. 

For  a  given  lift  coefficient  C^  passed  through  COMMON  USER  and  a  given 
Mach  number  M  transmitted  through  COMMON  B,  subroutine  GETALF  will  return  the 
corresponding  angle  of  attack  i . 

In  the  example  of  Figure  7.5  starting  with  an  initial  estimate  of  angle 
of  attack  a  Newton  Rahpson  method  is  used  to  evaluate  successive  approximations 
of  lift  coefficient  obtained  from  a  call  to  SURF  which  returns  ^(“.M). 
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7.2.6  AERO  Subroutine 


AERO  is  a  user  dependent  routine  called  by  TABLE  or  SEP  that  returns 
through  COMMON  B,  values  of  total  drag  and  effective  thrust. 

In  the  example  given  in  Figure  7.6  an  initial  estimate  of  the  angle  of 
attack  is  made  and  a  Newtonian  iteration  is  used  to  caluculate,  for  a  given  Mach 
number,  altitude  and  load  factor,  the  trimned  angle  of  attack  and  lift 
coefficient  using  the  curve  fitted  data  stored  in  array  E.  Effective  thrust  and 
total  drag  are  then  calculated. 

7.2.7  GETCD  Subroutine 

GETCD  is  a  routine  called  by  AERO  to  determine  the  drag  coefficient 
(including  incremental  store  drag)  for  a  given  lift  coefficient  and  Mach  number. 
This  can  be  done  in  various  ways  depending  on  the  data  available,  as  shown  m 
Figures  7.7  and  7.8. 

7.2.8  GETCDA  Subroutine 


The  purpose  of  GETCDA  is  to  return  the  drag  coefficient  and  angle  of 
attack,  given  Mach  number  and  lift  coefficient.  This  can  simply  be  done  bv 
calling  GETALF  and  GETCD  as  shown  in  Figure  7.9. 

This  routine  is  required  if  the  user  wishes  to  run  program  PC2P1  to  use 
the  spline  fitted  data  created  for  P1V2  to  make  up  an  aircraft  file  for  PACAM2. 

7.2.9  GETCDCL  Subroutine 


For  a  given  Mach  number  (passed  through  COMMON  B)  and  a  given  angle  of 
attack  (passed  through  COMMON  USER)  GETCDCL  determines  the  lift  coefficient  and 
then  the  drag  coefficient. 

In  the  example  shown  in  Figure  7.10,  Cp  is  determined  by  a  call  to  SURF 
to  access  the  first  spline  stored  in  array  E  representing  0^(1, M);  Cp  is 
then  evaluated  by  a  call  to  GETCD.  Another  example  as  shown  in  Figure  7.11. 

7.2.10  XSPMIN  and  XSPMAX  Functions 

XSPMIN  and  XSPMAX  are  used  in  all  the  user  defined  functions;  thev 
respectively  return  the  lower  and  upper  bound  of  X,  stored  in  the  two 
dimensional  spline.  These  functions  are  shown  in  Figure  7.12. 

7.2.11  YSPMIN  and  YSPMAX  Functions 


Equivalent  to  7.2.10  but  for  the  Y  variable. 


7.3  File  P1V2  (PI  Version  2) 


7.3.1  P1IN  Subroutine 

P1IN  controls  all  input  data.  All  ai rc raft - independent  parameters 
defining  the  type  of  calculations  required  are  accepted  in  conversational  mode, 
with  validity  checks  where  possible.  Aircraft-dependent  data  is  read  from  disk 
via  a  call  to  IDENT. 

This  subroutine  has  been  slightly  modified  to  allow  for  the  additional 
options:  turn  rate  variable  and  fourth  energy  parameter.  The  modifications  are 

labelled  K3,  K4 ,  K5  in  the  program  listing. 

7.3.2  BAD  I  NT  Subroutine 

BADINP  types  out  an  error  message  on  the  user's  terminal.  It  is  used  to 
check  input  validity  during  conversational  input  operations. 

This  subroutine  has  not  been  changed. 

7.3.3  TABLE  Subroutine 

TABLE  controls  the  calculation  and  output  of  the  unoptimised 
manoeuvrability  grid  using  the  following  grid  variables: 

.  '  being  energy  state  or  pressure  altitude  depending  on  the  input 

request . 

.  Mach  number. 

.  i  being  load  factor  or  turn  rate  depending  on  the  input  request. 

The  principal  outputs  are  now: 

.  Energy  parameter  :  Ps  or  Ps/wf  or  PgV/wf  or  Ps/V. 

.  Turn  rate  (if  -  is  load  factor)  or  load  factor  (if  C  is  turn  rate"!. 

.  The  load  factor  and  turn  rate  corresponding  to  the  maximum  lift 

boundary. 

The  subroutine  has  been  modified  in  Version  2  to  allow  for  the  additional 
options  (turn  rate  grid  variable  and  fourth  energy  parameter)  and  to  calculate 
specific  excess  power  Ps  using  gross  thrust  rather  than  net  thrust  (see 
Chapter  3  for  more  details).  The  modifications  are  labelled  K2,  K3,  K4,  K5  in 
the  program  listing. 

7.3.4  ATMOS  Subroutine 


ATMOS  calculates  atmospheric  parameters,  giving  temperature  as  a  function 
of  pressure  height  in  either  ICAO  or  ARDU  Atmospheres. 

This  subroutine  has  not  been  changed. 

7.3.5  INTRP  Subroutine 


INTRP  is  a  linear  interpolation,  returning  ordinate  and  gradient 
information. 

This  subroutine  has  not  bee  changed. 


*  i 
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7.3.6  HEIGHT  Subroutine 

HEIGHT  calculates  the  geopotential  height  corresponding  to  a  given 
pressure  height  in  the  nominated  atmosphere. 

This  subroutine  has  not  been  changed. 

7.3.7  HTRl'E  Subroutine 

HTRUE  calculates  the  true  height  for  given  atmospheric  parameters. 

This  subroutine  has  not  been  changed  and  is  explained  more  fully  in 
Reference  1 . 

7.3.8  IDENT  Subroutine 

IDENT  is  an  input  subroutine  which  reads  the  aircraft  identification 
data  file. 

It  has  been  modified  to  allow  the  option  of  either: 

.  calling  BINIh'  to  read  binary  records  from  an  engine  data  file  into  array 
F  and  an  aerodynamic  data  file  into  array  E. 

.  calling  PAC2IN  to  read  a  PACAM2  ASCII  data  file  (engine  and  aerodynamic 
tables!  into  arrays  E  and  F.  (See  Paragraphs  7.1  and  7.3  for  more 
details) 

7.3.9  P 1  OUT  Subroutine 

P10UT  provides  output  control  for  unoptimized  performance  calculations. 

A  text-formatted  output  file  is  produced  for  printing  and  an  alphanumeric  output 
file  is  produced  for  input  to  program  P2V2. 

This  subroutine  has  been  modified  to  print  the  lift  limit  boundaries,  and 
to  allow  for  the  extra  option  of  the  energy  parameter  Ps/V  and  the  alternative 
option  of  a  turn  rate  grid.  The  modifications  are  labelled  K2,  K3,  K4 ,  K5  in 
the  program  listing. 

7.3.10  BININ  Subroutine 

BININ  reads  data  into  a  nominated  storage  area  from  a  specified  disk 
file.  This  file  must  consist  of  a  sequence  of  binary  records,  with  the  first 
word  in  each  record  being  the  number  of  data  items  in  that  record. 

This  subroutine  has  not  been  changed. 

7.3.11  A1TIT  Subroutine 

ALTIT  performs  a  Newton  iteration  to  determine  pressure  height  hp  and 
other  atmospheric  quantities,  when  the  independent  variables  are  energy  state 
E6  and  Mach  number. 

This  subroutine  has  not  been  changed  and  is  described  in  detail  in 
Reference  1. 
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7.3.12  PARAMS  Subroutine 

PARAMS  calculates  airspeed  and  pressure  parameters. 

The  initial  angle  of  attack  estimate  for  each  aircraft  trim  calculation 
has  now  been  removed  from  this  subroutine  and  added  to  the  user  defined 
subroutine  AERO. 

7.3.13  MAXMAN  Subroutine 


MAXMAN  controls  the  calculation  and  output  of  the  optimised 
manoeuvrability  grid,  using  energy  height  Es  and  load  factor  or  turn  rate  as 
grid  variables.  Energy  parameter  (Ps,  Ps/wf,  PsV;wf  or  Ps^T^  ls 

optimised  by  varying  Mach  number,  and  a  simple  Euler  integration  is  performed  at 
each  energy  state  to  estimate  a  climb  time  historv. 

This  subroutine  has  been  modified  to  allow  for  the  additional  option  of 
the  energy  parameter  Ps/V  and  the  option  of  a  turn  rate  grid  variable.  The 
modifications  are  labelled  K5  and  K8  in  the  program  listing. 

7.3.14  SEP  Subrout  me 

SEP  calculates  the  energy  parameter  {P  ,  Pc/'w'f »  PgV/wf  or  ?s/7) 
for  the  optimised  grid. 

TK i s  subroutine  has  been  largely  rewritten.  The  calculation  of  Ps  has 
beer  improved  and  the  subroutine  now  allows  for  the  additional  options 
available.  The  mod i f icat ions  are  labelled  K3,  K3,  KS  in  the  listing. 

'.3.13  PiOVTA  Subroutine 


PlOl'TA  provides  output  control  for  optimized  performance  calculations,  in 
a  like  manner  to  subroutine  P10VT. 

This  subroutine  has  also  been  modified  to  allow  for  the  additional 
options  available.  The  modifications  are  labelled  K3,  K3,  K8  in  the  program 
list  i  r.  g  . 

7.3.16  MONSEP  Subroutine 


MONSEP  is  a  routine  called  by  ROMIN  to  monitor  the  convergence  of  the 
energy  parameter  optimisation. 

This  routine  has  not  been  altered  and  is  described  in  detail  in 
Reference  1 . 
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7.4  Subroutine  Library  P1LIB 


This  subroutine  library  is  exactly  the  same  as  described  in  Reference  1. 
It  contains  routines  for  curve  fitting  and  optimisation  which  are  explained  in 
detail  within  the  program  listing  and  are  described  very  briefly  here  for 
compl et  eness . 

7.4.1  Sl'RF  Subroutine 


SURF  takes  spline  fit  information  representing  a  surface  stored  in  vector 
form  and  provides  at  the  grid  point  (x,y)  the  function  z  =  f(x,y)  and,  if 
requested,  its  1st  and  2nd  order  derivatives.  The  spline  data  represent  knots 
and  B-splme  coefficents  for  least  square  curve  fits  of  z  against  x  for  discrete 
values  of  v . 

7.4.2  CUS1CS  Subroutine 


CL’BICS  applies  a  least-squares  cubic  spline  fit  of  B-sp’ines  to  weighted 
data  points  with  selected  knots.  In  the  present  task,  CUB1CS  is  only  required 
to  fit  a  simple  cubic  to  four  points;  hence  the  fit  is  exact,  all  weights  are 
unity,  and  no  external  knots  are  required. 

7.4.3  CRECKO  Subroutine 

CHECKD  checks  the  validity  of  data  supplied  to  subroutine  Cl'BICS  for 
curve  fitting. 

7.4.4  SPPER3  Subroutine 

SPDER3  evaluates  the  cubic  spline  f(x)  and  its  first  derivative  d  f  ( >; '  «'dx , 
based  on  normalised  B-spline  coefficients  and  associated  knot  positions. 

7.4.5  ROM!  N  Sub  rout  me 

ROMIN  finds  the  local  minimum  of  an  unconstrained  function  of  n 
variables. 

7.5  Subroutine  Librarv  PAC2IN 


This  subroutine  librarv  contains  utility  routines  for  use  when  the  data 
files  provided  are  in  ASCII  tabular  format. 

7.5.1  PAC2IN  Subroutine 

PAC2IN  is  a  new  routine  which  has  been  added  to  read  the  relevant 
information  of  a  PACAM2  aircraft  data  file  which  consists  basically  of  valuer  of 
military  and  after-burner  thrust,  fuel  flow  and  ram  drag  stored  in  array  form 
with  Mach  number  and  altitude  as  arguments.  Aerodynamic  data  in  the  form  of 
drag  coefficients  and  angle  of  attack  are  also  read  in  using  Mach  number  and 
lift  coefficient  as  arguments. 
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Hence  PAC2IN  first  reads  in  the  number  of  altitude,  Mach  number  and  lift 
coefficient  arguments  for  which  data  will  be  read.  The  actual  values  of  the 
altitudes,  Mach  number  and  lift  coefficient  arguments  are  then  read  followed  by 
the  maximum  Mach,  number  for  each  altitude  argument  and  the  maximum  for  each 
Mach  number  argument.  Tabulated  values  of  military  thrust,  after-burner  thrust, 
military  fuel  flow,  after-burner  fuel  flow,  and  intake  momentum  drag  are  then 
read  for  the  values  of  the  altitude  and  Mach  number  arguments  listed  and  entered 
into  array  F  via  calls  to  TABIN. 

Finally  values  of  drag  coefficient  and  angle  of  attack  for  the  listed 
values  of  the  M  and  Cj^  arguments  are  entered  into  array  E  via  calls  to  TABIN. 

7.5.2  TABIN  Subroutine 

TABIN  transfers  data  into  a  nominated  internal  storage  area  from  a 
specified  disk  file.  Although  presently  this  file  is  in  ASCII  and  has  the 
format  of  a  PACAM  aircraft  data  file  (doubly  dimensioned  table),  the  adaptation 
to  a  different  tabulated  format  could  be  easily  accomplished. 

7.5.3  PACINT  Subroutine 

PACINT  is  a  double  linear  interpolation  routine.  Given  a  variable  (e.g. 
thrust)  that  is  tabulated  as  a  function  of  two  independent  arguments  (e.g. 

Mach  number  and  altitude),  a  value  for  this  variable  at  any  set  of  entries  is 
determined  by  double  linear  interpol at  ion  between  adjacent  tabulated  points. 

7. 5. A  FIND  Subroutine 

FIND  will  return  the  position,  within  a  single  dimensioned  ordered  arrav, 
of  the  tabular  argument  above  the  input  variable. 

7.5.5  DINTRP  Subroutine 

Given  z,.  =  flXj,  y^)  zJ0  =  f(Xj,  y-J 

z , ,  =  f  ( x ,  y  )  z-,0  =  f  (x.,,  y.,) 


DINTRP  will  perform  a  double  linear  interpolation  and  return  z  =  f(x,v). 
7.5.6  PACDER  Subroutine 

For  a  variable  z  that  is  tabulated  as  a  function  of  two  independent 
arguments  x  and  v,  PACDER  determines  a  value  for  z  and  its  derivatives  with 
respect  to  x  and  v  a:  any  set  of  entries  (x,y). 


SUBROUTINE  THRUST 


C  WRITTEN  BY  VR  ARL  1983. 

C  THRUST  RETURNS  FOR  CIVEN  HFIOHT  /MACH  NO.  - 
C  (1)  WF  -FUEL  CONSUMP.  (  TO  PI  THRU  COMMON  fl /  3 

C  (23  XMLIM-HAX.  MACH  ALL.  (  .  "  "  ) 

C  (3)  FN  -NET  THRUST  (  ONLY  TO  AFRO.AEROMX  THRU  COMMON  USF'M 

C  (4)  FC  -CROSS  THRUST  (  . 

C  (5)  RD  -RAM  DRaC  (  . 

C  MUST  BE  CALLED  REFORE  EITHER  AERO  OR  AEROMX  ARE  CALLED 


PARAMETER  NUMENC- 1 
COMMON  /  B/  BCOni 
COMMON  ,'C!  C(200) 

COMMON  / FDATA/  F(2000) 

COMMON  /USER/  FN  ,FO,RD,CLT  , CD  .ALPHA .  CDS 


EQUIVALENCE 

( 

XM, 

i  3( 

33). ( 

HP , 3C  o )  3 , C 

7 ,  R  (  Q  ) 

1 

( 

PLA, 

,B( 

153), 

1 

( 

SPEF , 

,B( 

183). ( 

WF,B(  42)). ( 

3 

< 

XMLIM, 

,B( 

61)3. ( 

IERR , R(  75)) 

EQUIVALENCE  (  C,C( 39 ) ) 


C  ORDER  OF  SPLINES  IN  F  IS  ASSUMED  TO  RE  - 
C  (1)  NET  THRUST-It.T-  (MACH  .HEIGHT )  I LRF  ] 

C  (23  ENCINE  RAM  DRAC-INT-  (MACH  .HEIGHT  3  (LRF’ 

C  (3)  FUEL  CONS.-INT-  (MACH , HEIGHT  3  ! LRM/HR 1 

C  (43  NET  THRUST -A/ 8-  ( MACH .HEIGHT 3  1 LRF ] 

C  (5)  ENCINE  RAM  DRAG-A/R-  (MACH  .  HF.ICHT  3  I  LRF  1 

C  '  (6)  FUEL  CONS. -A/ B-  (MACH, HEIGHT  3  | LRM/HR 1 


C  CHECK  IF  A/B  OR  INTERMEDIATE  POWER  REOtTSI  - 
C  IF  NOT  STOP  WITH  APPROPRIATE  MESSAGE. 

NJ  *  2 

IF(  AB  S(  PLA-200 .3.LE..0O13  GOTO  10 
IF( ABS( PLA-100. ) .LE. .001 3  GOTO  20 
TYPE  25  . FLA 

75  FORMAT  (  '  CANNOT  PROCESS  '.F5.2,'  T  THRUST  REOUFST  YET '  3 

STOP 

C  A/B  POWER  SO  READ  THRU  TO  START  OF  A/B  SPLINTS 

10  DO  11  1-1,3 

NJ  -  NJ  +  F(  NJ-1  )  *  .00] 

11  CONTINUE 

C  IF  INTERMEDIATE  POWER  REOUEST  ALREADY  AT  START  OF  INT  SPLINES 


C  SET  NOMINAL  MACH  NO.  LIMIT 
20  XMLIM  -  XSPMAXC  F(NJ)  ) 


C  FIND  NET  THRUST.FN. 

CALL  SURF(F(NJ),XM.HP,FN,D1,D2,D3,  IERR  ) 
IF(  IERR.  NF..03  GOTO  900 
FN  -  FN»NUMENG 

C  FIND  RAM  DRAG ,  RD . 

NJ  -  NJ  +  F(NJ-l)  +  .001 

CALL  SURF(  F(NJ)  .XM.HP.RD.Dl  .D2.D3,  IFRR3 

IF( IERR. NE.03  COTO  900 

RD  ■  RD*NUMFNC 

FC  -  FN  +  RD 

C  FIND  FUEL  CONSUMPTION, WE. 

NJ  -  NJ+F ( NJ-I 3+.001 

CALL  SUR  F(F(NJ3,XM,HP,WF,D1,D2,D3,  IER  R  3 
IF ( IERR. NE.03  COTO  900 
WF  -  WF  *  NUMENC 

900  CONTINUE 
RETURN 
END 


FIG.  7.1  Example  1  of  Subroutine  THRUST 
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SUBROUTINE  THRUST 
C  WRITTEN  BY  VR  ARL  MARCH  1983 
C  THRUST  RETURNS  FOR  CIVEN  KFICHT  /MACH  NO.  - 
C  (1)  WE  -FUEL  CONSUME.  (  TO  PI  THRU  COMMON  /R/  ) 

C  (2)  XMLIM-HAX.  MACH  ALL.(  .  "  "  ) 

C  (3)  FN  -NET  THRUST  (  ONLY  TO  AERO.AFROMX  THRU  COMMON  /USFR/) 

C  (4)  FC  -CROSS  THRUST  ( .  "  "  “  ) 

C  (5)  DELTL-LIFT  LOSS  (  “  "  "  M  "  "  "  1 

C  (63  BD  -RAM  DRAC  ( .  "  "  "  ) 

C  MUST  BE  CALLED  BEFORE  EITHER  AERO  OR  AEROMX  ARF.  CALLED 

PARAMETER  AJ  -  6.38  '.NOZZLE  AREA  (FT**  2) 

REAL  JVR. LRR.K1 ,  K2 
COMMON  / B/  B( 200 ) 

COMMON  /FDATA/  F(2000) 

COMMON  /USER/  FN. FC.RD.CLT, CD, ALPHA, DELTL .THETAE , 

1  CLO,  DCLDA,  K1 ,  K2  ,CD0,  CLCRIT  .CDS 

EQUIVALENCE  (  XM,B(  3)).(  T,B(  4)),(  HP,B(  6)), 

1  (P  ,8(  7)3, (  AJN,R(  13)), 

1  (  SREF , B(  16)), (  UF,B(  42)), (  DIM,R(49  )), 

3  (  XM LIM , B (  61)), (  IERR , R(  75)) 

C  ORDER  OF  SPLINES  IN  F  IS  ASSUMED  TO  BE  - 
C  Cl)  (  NORMALIZED  ]  CROSS  THRUST -COMRAT-  <  MACH  .HF1CHT  ) 

C  (2)  [ NORMALIZED ]  FUEL  CONS. -COMBAT-  (MACH .HFICHT ) 

C  (3)  [NORMALIZED]  NET  THRUST -COMRAT-  (  MACH  ,  HFICHT  ) 

C  (4)  LIFTLOSS/FG  RATIO  (JVR.ANC  OF  NOZZLE  DFF ) 

C  POWER  LEVEL  FOR  AIRCRAFT  IS  ASSUMED  TO  ALWAYS  BE  COMBAT. 

C  CALCULATE  NORMALIZING  FACTORS  -  PDPO  -  p/p 0,  TtTTO  -  T/TO 
PDPO  -  P/2116.22 
TDTO  -  T/288. 15 

C  SET  NOMINAL  MACH  NO.  LIMIT 

XNLIM  -  XSPMAXt  F(2>  ) 

C  FIND  CROSS  THRUST  FC  FROM  SPLINE  THEN  UN-NORMALI SF. 

NJ  -  2 

CALL  SURFCF(NJ),XM,HP.FG.Dt.D2,D3,IERR) 

IF(IERR.NE.O)  GOTO  900 
FC  -  FC*PDPO 

C  FIND  FUEL  CONSUPTION  WF  FROM  STLINE  THEN  UN-NORMALI  SF 
NJ  -  NJ*-F(SJ-1  H.OOl 

CALL  S(/RF(F(NJ),XM,KP,WF,  31.02,03,  TERR) 

IF ( IERR. ME. 0)  GOTO  900 
WF  -  WF*PDP0*S0RT(TI7r0) 

C  FIND  NET  THRUST  FROM  SPLINT.  THEN  UN-NORMALI  SEFN . 

NJ  -  NJ+F(NJ-1  J+.OOl 

CALL  SURF(F(NJ),XM.HP,FN,D1,D2,D3,  IERR) 

IF(  IERR..VE.O)  GOTO  900 
FN  -  FN*PDPO 

C  CALCULATE  RAM  DRAC  RD 
RD  -  FC  -  FN 

C  FIND  LIFT  LOSS  RATIO(LRR)  4  CALCULATE  LIFT  LOSS(DELTL) 

DELTL  -  0.0 

IF  (  XM .  LT .  IE-4)  GOTO  900 
NJ  -  NJ  +  F(NJ-l)  +  .001 

JVR  .  SORT(( FC* SREF)/ (2.0*DIM*AJ ) )  'JET  VEL.  RATIO 

JVR  -  AMINKJVR, XSPMAX(  F(N.I)  )) 

AJN 1 -AM  IN  1  ( A  JN ,  YSPMAX(  F(  NJ  )  )  ) 

CALL  SURF(F(NJ),JVR,AJN1,LRR,D1,D2,D3,  IERR) 

IF(IERR.NE.O)  GOTO  900 
DELTL  -  LRR  *  FC 
900  CONTINUE 
RETURN 
END 


FIG.  7.2  Example  2  of  Subroutine  THRUST 
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SUBROUTINE  THRUST 

C  THRUST  RETURNS  FOR  GIVEN  HE  I  HIT  /MACH  NO.  - 
C  (i;  WF  -FUEL  CONSUMP.  <  TO  PI  THRU  COMMON  / H/  ) 

C  (2?  XMLIM-MAX.  MACH  ALL.  (  .  "  "  ) 

C  (31  FN  -NET  THRUST  (  ONLY  TO  AERO  .AEROMX  THRU  COMMON  /UEF.F/) 

C  Hi  FC  -CROSS  THRUST  ( . .  "  "  "  ) 

C  (51  RD  -RAM  DRAC 

C  (6)  IXM  -INDEX  S.T.  VMARC(IXM-l)  <  XM  <«  VMARC(IXM) 

C  MUST  BE  CALLED  BEFORE  EITHER  AERO  OR  AFROMX  ARE  CALLED 

COMMON  /B/  8(200) 

COMMON  / FDATA/  FC2000) 

COMMON  /USER/  U(10) 

DIMENSION  HARCCO),  VHARCCO ),  VMTOPCO) 

EQUIVALENCE  (  XM,B<  3)),(  T,R(  4)).(  HP,B(  A)), 

1  (P  ,B(  7)),<  AJN.Bf  13)). (  PLA ,  B  (  15)), 

1  (  SREF.Bf  16)), (  WE,R(  42)1,  (  DIM,R<40  )), 

3  (  XMLIM,8(  61)), (  IERR, B(  75)) 

EQUIVALENCE  (  FN , U(  1)).(  FC.U(  2) ) ,(  RD.IK  3)), 

1  (  IXM.U(IO)) 

EQUIVALENCE  (  NA.FO  )),(  NM,F(  2)),(  HARC,F(  3)), 

1  (  VMARC,F(23)),(  VMTOP , F( 43) ) 

FINTRP(X.X1,Y1 ,X2,Y2)«YI+(X-X1 )*(  Y2-Y1 )/(X2-X!) 

;  assumes  a  paca.m2  a/c  data  file. 

;  assumes  f  contains  - 

:  (1)  NA  -NO.  OF  ALTITUDE  ARCS 

C  (2)  MM  -NO.  OF  MACH  ARCS 

C  (3)  RARC  -ARRAY  OF  ALTITUDE  ARCS 

C  (4'  VMARC  -ARRAY  OF  MACH  ARCS 

C  (5)  VMTOP  -ARRAY  OF  MAX  MACH (  ALT) 

C  (6)  THRMP  -TABLE  OF  MIL.  THRUST  (MACH  .ALT  1 

C  (71  THRA8  -TABLE  OF  A/B  THRUST(  MACH  .ALT ) 

C  (?)  DINM  -TABLE  OF  RAM  DRAC(  MACH  ,  ALT  1 

C  (91  FCMP  -TABLE  OF  MIL.  FUEL  CONS.  ( MACH  .ALT'- 

C  (lO)FCAR  -TABLE  OF  A/B  FUEL  CONS.  (  MACH  , ALT) 

C  FIN'D  IXM.  IHP 

CALL  FIN'D!  'MARC. NM.XM,  IXM  1 
CALL  F  IND(  HARG ,  NA,  HP  ,  IHP  1 

C  FIN'D  MAX  MACH  NO.  FOR  THIS  ALTITUDE 

XM. LIM  -  FINTRP(HP.KARC(  IHP-l  )  ,7MT0P(  IHP-1  ),HARC(  IMP)  , 

1  VMTOP(IHP)) 

C  FIN'D  CROSS  THRUST  :FC 

NJ  -  63  !  POINTER  TO  CURRENT  TABLE  IN  F 

FCMIL  ■  PACIMT(XM, HP. VMARC, HARC.F(SJ),  IXM,  IHP, NM.NA) 

NJ  •  NJ  *  NA«NM 
IFfPU.CT.lOO.)  GOTO  10 
C  INTERPOLATE  FROM  0  TO  MIL  POWER  RATING 
FG  -  FINTRP(PLA,0.,0., 100., FCMIL) 

GOTO  20 

C  INTERPOLATE  BETWEEN  MIL  POWER  6  A/B  POWER 

10  FGAB  -  PACINT(XM, HP, VMARC, HARC,F(NJ), IXM, WP, NM.NA) 

FC  -  FINTRP(PLA, 100. .FCMIL, 200., FGAB) 

C  FIND  RAM  DRAC  -  INTAKE  MOMENTUM  DRAC 
20  NJ  *  NJ  *  NA*NM 

RD  -  PACINT(XM,  HP,  VMARC,  HARC.F(NJ),  IXM,  IHP,  NM.NA) 

C  CALCULATE  NET  THRUST  :FN 
FN  -  FC  -  RD 

C  FIND  FUEL  CONSUMPTION  :WF 
NJ  -  NJ  +  NA*NM 

WFMIL  •  PACINTtXM,  HP,  VMARC,  HARC.F(NJ),  IXM,  DIP,  NM.NA) 

IF( PLA. CT . 100. )  COTO  50 
C  INTERPOLATE  BETWEEN  0  6  MIL  POWER 

WF  -  FINTRP(PLA,0. ,0. , 100. , WFMIL) 

OOTO  100 

C  INTERPOLATE  BETWEEN  MIL  POWER  6  A/ B  POWER 
50  NJ  -  NJ  +  NA*NM 

WTAB  -  PACINT(XM,  HP,  VMARC,  HARG,F(NJ),  IXM,  IHP,  NM.NA) 

WF  •  FINTRP(PLA, 100. .WFMIL, 200., WFAB) 

C 

100  RETURN 

END 


FIG.  7.3  Example  3  of  Subroutine  THRUST 
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SUBROUTINE  AEROMX 


C 

C  AEROMX  FINDS  FOR  A  CIVEN  HEIGHT  &  MACH  NO. 

C  (1)  CLMAX 

C  (2)  ALPMAX  -  ALPHA  CORRESPONDING  TO  CLMAX 

C  (3)  TLMAX  -TOTAL  LIFT  CORRESPONDING  TO  CLMAX 

C  (4)  CDS  -STORE  DRAC  COEF.<  PASSED  THRU  /USER/  TO  AFRO) 

C  AEROMX  ASSUMES  A  CALL  HAS  BEEN  MADE  TO  THRUST  PRIOR  TO  ITS  CALLING. 

C  AEROMX  MUST  BE  CALLED  BEFORE  AERO  IS  CALLED  FOR  THAT  PARTICULAR  HEIGHT/ 

C  MACH  NO.  COMBINATION. 

C  VR  ARL  APRIL  1983 

C  ASSUMES  ORDER  OF  SPLINE  FITS  IN  F  IS  - 
C  (1)  CL< ALPHA, MACH) 

C  (2)  CLMAX(MACH) 

C  (3)  CDL(CL.MACH) 

C  (4)  CDO(MACH) 

COMMON  /B/  B( 200) 

COMMON  /C/  C(200) 

COMMON  /EDATA/  EC2000) 

COMMON  /USER/  FN ,  FC.RD.CLT , CD, ALPHA,  CDO, CDS 
DIMENSION  XMDCDS(IO),  DELCDS(20) 

EQUIVALENCE  ( RADIAN, C(  40)) 

EQUIVALENCE  (  GN,B(  2)),(  XM,B(  3)),(  THETAI,R(  13)), 

1  (  W,B<  17)), (  EFFT , B( 4 1 ) ) , 

1  (  D,B(  43)), (  CLMAX, B(  44)), 

2  (ALPMAX, B(  45)), (  TLMAX, B(  48)), (  DIM,B(  49)), 

3(  IERR, B(  75)). (  NDCDS,B(  80) ) , ( XMDCDS, B(  81 ) ) , ( DELCDS, B( 1 0 1 ) ) 

FINTRP(X,X1,Y1,X2,Y2)  -  Y1  +  (X-Xl)  *  (Y2-Y1)  /  (X2-X1) 

C  FIND  CLMAX  FROM  2ND  SPLINE  FIT  IN  E 

IERR  •  0  !  RESET  IERR 

NJ  -  E(  I  H2.0001  !NJ  POINTS  70  2ND  SPLINF  IN  E 

X  -  AMAX1 (XM,XSPMIN(  E(NJ)  )) 

Y  -  0 

CALL  SURF(E(NJ), X,Y, CLMAX,  D1,D2,D3,  IERR)  !  CET  CLMAX  EOF  MACH  X 

IF(IERR.NE.O)  GOTO  30 

C  FIND  CDO  FROM  4TH  SPLINT.  IN  E 
NJ-NJ+E(NJ-1  )+.0Ol 
NJ-NJ+E(NJ-1  )+.00l 
X-AMAX1 ( XM, XSPMINC  ECNJ ) ) ) 

Y-0 

CALL  SURF(E(NJ) ,X, Y.CDO, Dl , D2, D3, IERR) 

C  FIND  ALPHA  CORRESPONDING  TO  CLT  -  CLMAX 
CL*  -  CLMAX 
CALL  GETALF 
ALPMAX  -  ALPHA 

C  CALCULATE  TLMAX  :  LIFT  CORRESPONDING  TO  CLMAX 
TLMAX  -  FC*  SIND  (ALPMAX)  +•  CLMA.X*DI.M 
30  Cnf.TINUE 

C  FIND  STORE  DRAC 
CDS  •  0. 

IF  (NDCDS.EO.O)  GOTO  900 
IF  (XN.LE.XMDCDS(l))  CDS  -  DF.LCDS(l) 

IF  (XM.CD.XMXDS(NDCDS))  CDS  -  DELCDS(NTXIDS) 

IF  (XM.LE,XMDCDS(  1  )  .OR.  XM.  GE .  XMDCDS  ( NDCDS  ) )  GOTO  540 
DO  450  1*2 ,  NDCDS 

450  IF  (XM.LE.XMDCDS(I))  GOTO  500 
IERR  -  104 
GOTO  900 
500  Dt  -  1-1 

CDS  -  FINTRP  (XM,  XMDCDSUM),  DELCDS(IM),  XMDCDS(I),  DF.LCDS(I)) 

540  CDS  -  CDS  *  1.0E-04 
900  CONTINUE 
RETURN 
END 


1  NJ  POINTS  TO  4TH  SPLINF  IN  F 

! GET  CDO  FOR  MACH  X 


FIG.  7.4  Example  of  Subroutine  AEROMX 


1 


SUBROUTINE  GETALF 

C  FINDS  ALPHA  FOR  GIVEN  CLT  IN  /USER/  AND  XM  IN  /B/ 

C 

COMMON  /EDATA/  E(2000) 

COMMON  /B/  B( 200 ) 

COMMON  /USER/  FN ,FG,RD , CLT , CD .ALPHA , CDO, CDS 
EQUIVALENCE  (XM  ,3(  3)!,(  IERR,B(  75)) 


C  USE  N/R  TO  FIND  FROM  CL (ALPHA, MACH)  -THE  1ST  SPLINE  FIT  IN  E  - 
C  THE  ALPHA  CORRESPONDING  TO  CLT. 

ALPHA  =  CLT  /.09  l£WP.  DET  BEST  1ST  ESTIMATE 

NJ  =  2  !NJ  POINTS  TO  1ST  SPLINE  IN  E 

ASPMAX  =  XSPMAX( E(NJ) ) 

YM  =  AMAX1 (XM, YSPMINC  E(NJ)  )) 

DO  20  1=1,20 


IERR  =  -1  IDERIVATIVES  REQUIRED 

A  =  AMINKASPMAX, ALPHA) 

CALL  SURF( E(NJ ) ,A , YM,CL, DCDA, D2, D3, IERR ) 

IF(IERR.NE.O)  GOTO  21 

CL  =  CL«(ALPHA-A)»DCDA  ! EXTRAPOLATE  OUTSIDE  RANGE 

FA  =  CL-CLT 

IF(ABS(FA).LT.1E-V)  GOTO  25  IEXI?  AS  CLOSE  ENOUGH 
ALPHA  =  ALPHA  -  FA/DCDA 


20 

CONTINUE 

21 

IERR  =  102 

I  DID  NOT  CONVERGE  IN  20  ITER 

25 

RETURN 

END 

FIG.  7.5 


Example  of  Subroutine  GETALF 
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o  o  o  o  o  o  o 


SUBROUTINE  AERO 


AERO  ROUTINE  RETURNS  FOR  CIVEN  LOAD  FACTOR  /HEIGHT  /MACH  NO.  - 

(1)  D  -  DRAG 

(2 )  EFFT  -  EFFECTIVE  THRUST 

AERO  ASSUMES  AEROMX  HAS  PREVIOUSLY  BEEN  CALLED  TO  SET  UP  VARIABLES  DEPENDING  ONLY  ON 
GIVEN  HEICHT/MACH  NO. 

SIMILAIRLY  IT  ASSUMES  THRUST  HAS  SET  UP  PROPULSION  ASSOCIATTED  VAR  FOR  THE 
GIVEN  HEIGHT  /MACH  NO. 


C  ASSUMES  ORDER  OF  SPLINE  FITS  IN  E  IS  - 
C  (1)  CL (ALPHA, MACH) 


COMMON  /B/  8(200) 

COMMON  10/  C(200) 

COMMON  /EDATA/  E(2000) 

COMMON  /USER/  FN ,FG,RD, CLT, CD, ALPHA, CDO, CDS 
DIMENSION  XMDCDS(20),  DELCDS(20) 

EQUIVALENCE  ( RADIAN ,C(  VO)) 

(  GN,B(  2) ) . (  XM,B(  3)), 
(  H,B<  17)). (  EFFT,B(  VI)), 


EQUIVALENCE 

1 

1 

2 

3 


(THETAI,B(  13)), 


D,B(U3  )),(  CLMAX, B(  uu ) > , 
(LMAX.Bl  V8)),(  DIM,3(  V?)),( 


lERB.Bt  75)), 


(  NDCDS ,B(  80)) ,(XMDCDS,3(  31 ) ) , (DELCDS ,B( 101 ) ) 


FINTRP(X,X1,Y1 ,X2,Y2)  :  II  >  (X-X1 )  •  (Y2-Y1)  /  (X2-X1  ) 


C  FIND  TRIMMED  ANCLE  OF  ATTACK  A  CL  FOR  GIVEN  HIEGHT ,MACH  A  LOAD  FACTOR. 


C  FIND  CL  TRIM  (CLT)  A  ALPHA  USINC  N/R  .ITERATING  ON  ALPHA 

T1  =  GN»W  ’TEMP.  VAR  TO  SAVE  CALC. 

IFtXM.GE.O.O)  GOTO  110 
IERR  =  103 
GOTO  900 


110 


NJ  =  2  ! POINTS  TO  1ST  SPLINE 

DCLDA  =  0.00  ! ESTIMATE  FROM  GRAPH 

ALPHA  «  T1  /  C  FG/RADIAN  *  DIM*DCLDA  )  ! INITIAL  EST. 

ASPMAX  =  XSPMAXt  E(NJ)  ) 

DO  120  1=1,10 

IERR  =  -1  ! REQUIRE  SPLINE  DERIV3 

A  =  AMINKASPMAX,  ALPHA) 

CALL  SURF(E(NJ) , A, YM, CLT, DCLDA, D2.D3, IERR) 

IF(IERR.NE.O)  GOTO  9 00 

CLT  =  CLT»( ALPHA-A )*DCLDA  I  LINEAR  EXTRAPOLATION 

DFDA  =  DIM'DCLDA  *  FC»COSD( ALPHA ) /RADIAN 
F  =  FG*SIND< ALPHA)  »  CLT*DIM  -  T1 
DALF  =  F /DFDA 

IF(ABS(DALF).LT. 0.001)  GOTO  200  ! EXIT  AS  CLOSE  ENOUCH 

ALPHA  =  ALPHA  -  DALF 


120  CONTINUE 

IERR  =  102  IMORE  THAN  20  ITERATIONS 

GOTO  900 


C  CALCULATE  EFFT  -  EFFECTIVE  THRUST 
200  EFFT  =  FG*COSD( ALPHA )-RD 


C  CALCULATE  CD  4  THEN  D 
CALL  GETCD 
D  =  CD  •  DIM 
900  CONTINUE 

RETURN  I  RETURN  NORMALLY 

END 


FIG.  7.6  Example  of  Subroutine  AERO. 
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SUBROUTINE  GETCD 


C  ROUTINE  TO  CALCULATE  CD  GIVEN  CLT  A  XM  AND  A  PRIOR  CALL  TO  AEROMX. 

COMMON  /USER/  FN ,FC, RD , CLT , CD .ALPHA , DELTL .THETAE, 

1  CL0.DCLDA.K1 , K2 .CDO.CLCRIT .CDS 

REAL  K1.K2 

C  CALCULATE  DRAG  COEFF 

IFtCLT.LE.CLCRIT!  CDL  s  (CDO.X1«CLT*CLT) 

IF(CLT.GT.CLCRIT)  CDL  »  ( CDO+(K 1-K2)*CLCRIT*CLCRIT 
1  -K2*CLT»CLT) 

C  CALCULATE  TOTAL  DRAG 
CD  *  CDL  CDS 
CONTINUE 
RETURN 
END 


Example 


of  Subroutine  GETCD 


SUBROUTINE  GETCL 

C  CALCULATES  CD  ' =  CDL<-CDO-CDS '  FnR  GIVEN  CLT  4  MACH. 

C  ASSUMES  CDL, CL. MACH  IS  THE  >RD  SPLINE  IN  E 

COMMON  / EDATA/  EI2000 
COMMON  .'B/  31200' 

COWON  USER  FN.FG.RD.  CLT.  CD.  ALPHA,  COO,  CDS 
EQUIVALENCE  XM  Bi  3)i..  IERR.Bi  ’S)> 

NJ  =  Ell.  <•  2.00) 

NJ  j  NJ  .  E(NJ-I)  .  .001  IPTS  TO  3RD  SPLINE  IN  E 

XCLT  j  AMINHCLT.XSPMAXIElN.,-'" 

YM  i  AMAX1  (  XM  ,  YSPMINt  E!NCi  )'. 

IERR  »  -1  f DERIVATIVE  REQUIRED 

CALL  SURF ( E!  NJ  ;  , XCLT . TM ,  CDL . DCDDCL . D2 , D3 , IERR ) 

IF(IERR.NE.O)  GOTO  TOO 

CDL  «  CDL*( CI.T-XCLT I •  DCDDCL  (LINEAR  EXTRAPOLATION 

CD  s  CDO  *  CDL  ♦  CDS 

900  CONTINUE 
RETURN 
END 


FIG.  7.8  Example  2  of  Subroutine  GETCD 
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SUBROUTINE  GETCDA 


C 

C  RETURNS  CD  4  ALPHA  FOR  OIVEN  XM  4  CL. 
C 

CALL  GETALF 
CALL  GETCD 
RETURN 
END 


FIG.  7.9  Example  of  Subroutine  GETCDA 


SUBROUTINE  GETCDCL 

C  RETURNS  CD  4  CL  FOR  GIVEN  MACH  4  ALPHA. 

COMMON  /EDATA/  2(2000) 

COMMON  /B/  B(200) 

COMMON  /USER/  FN ,F0,RD ,CLT .CD .ALPHA . DELTL.THETAE .CDDIFF , CDS 
EQUIVALENCE  (  XM,B(  3)M  IERR , 8(  T«)) 

NJ  =  2 

YM  =  AMAX1 (XM. YSPMINt  E(NJ))  ) 

CALL  SURF(£(NJ) .ALPHA, YM, CLT, D1.D2.D3. IERR) 

IF(IERR.NE.O)  RETURN 

CALL  GETCD 

RETURN 

END 


FIG.  7.10  Example  1  of  Subroutine  GETCDCL 


SUBROUTINE  GETCDCL 

C  RETURNS  CLT  4  CD  FOR  GIVEN  VALUES  OF  MACH  4  ALPHA. 

C  ASSUMES  PRIOR  CALL  TO  AEROMX. 

COMMON  /USER/  FN ,FG ,RD , CLT .CD, ALPHA . DELTL.THETAE , 
1  CL0,DCLDA,K1,K2,CD0,CLCRIT,CDS 

REAL  K1.K2 

C  CALCULATE  TRIWED  CL 

CLT  s  DCLDA*ALPKA  *  CLO 
C  CALCULATE  CD 

CALL  GETCD 

RETURN 

END 


FIG.  7.11  Example  2  of  Subroutine  GETCDCL 
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FUNCTION  XSPMIN(SPL) 


C 

C  RETURNS  LOWER  BOUND  OF  SPLINE  FIT  SPL  'S  I  RANGE 
DIMENSION  SPUD 

NC  *  SPL(2>  ♦  .001 
XSPMIN  =  SPL(NC»3) 

RETURN 

END 

FUNCTION  XSPMAX(SPL) 

C  RETURNS  UPPER  BOUND  OF  SPLINE  FIT  SPL’S  X  RANGE. 
DIMENSION  SPL  CD 

NC  r  SPL(2)  *  .001 
XSPMAX  =  SPL(NCUI) 

RETURN 

END 

FUNCTION  YSPMIN(SPL) 

C 

C  RETURNS  LOWER  BOUND  OF  SPLINE  FIT  SPL'S  Y  RANGE. 
DIMENSION  SPUD 

YSPMIN  -  SPL (3) 

RETURN 

END 

FUNCTION  YSPMAX(SPL) 

C 

C  RETURNS  UPPER  BOUND  OF  SPLINE  FIT  SPL'S  Y  RANGE. 
DIMENSION  SPUD 

NC  =  3PU2)  *  .001 
YSPMAX  =  SPL(NC*2> 

RETURN 

END 


FIG.  7.12 


Useful  Functions  used  in  ACRAFT 
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8. 


PROGRAM  P2V2  DESCRIPTION 


8 . 1  Loading  of  Program  on  DEC  System-10 

The  commands  for  loading  program  P2V2  and  saving  a  core  image  on  disk  on 
DEC  system-10  are: 

.  LOAD  P2V2,  P24LIB/ SEARCH 
.  SAVE  P2V2 

The  routines  contained  in  each  of  these  files  are  listed  in  the 
following  table  and  described  briefly  in  this  chapter. 


P2V2 

MAIN 

P2 IN  IN HMD  PSCON  PDSIFF 

RATE  1 

SEPCON 

RATE2 

MMP  GRID 

P24LIB 

INTRP 

BADINP 

INLAB 

PLTLAB  UNITS  WRLAB 

8.2  File  P2V2  (P2  Version  2) 


Most  subroutines  contained  in  this  file  have  been  extensively  rewritten. 
Their  purpose  is  described  here  and  further  details  are  available  in  the  source 
program  listing  which  is  extensively  documented. 

8.2.1  P2  MAIN  Program 

P2V2  processes  data  files  prepared  for  a  particular  aircraft  by  program 
ACRAFT  and  produces  either  plotter  files  ready  for  submission  to  the  off  line 
plot  queue  or  output  files  to  be  further  processed  by  the  contour  plotting 
program  P4V2 . 

Plotter  files  (named  P2.PLT)  are  produced  by  options  4B,  4C  and  4Z. 

The  output  file  P24A1 .CON  is  produced  by  option  4A  and  the  output  files 
P24A1.CON,  P24A2.C0N  and  P24DIF.CON  are  created  by  option  4D.  Option  4G  creates 
the  file  P24G.CON.  All  these  output  files  (.CON)  are  used  as  input  files  to  the 
contour  plotting  program  P4V2 . 

The  mam  program  has  been  modified  mainly  to  cater  for  the  additional 
option  4G.  The  Modifications  are  labelled  K1  and  K4  in  the  source  program 
1 ist ing. 

8.2.2  P2IN  Subroutine 


P2IN  is  called  by  the  mam  program  to  process  the  aircraft  identification 
header  on  the  input  file  and  requests  an  option  code  for  the  type  of  run.  An 
error  halt  is  forced  if  the  8000  words  allocated  for  vector  WORK  are 
insufficient  to  accommodate  the  data  required  by  the  input  grid.  It  is  the 
user's  responsibility  to  ensure  that  the  option  requested  is  consistent  with  the 
purpose  for  which  program  ACRAFT  produced  the  data. 

This  subroutine  has  been  modified  to  allow  for  the  additional  option  4C 
which  plots  contours  of  the  energy  parameter  on  a  Mach  number,  turn  rate  (or 
load  factor)  grid. 


8.2.3  INMMD  Subroutine 


Subroutine  INMMD  is  called  by  subroutine  P2IN.  It  performs  a 
conversational  dialogue  with  the  user  to  accept  text,  grid  data  and 
configuration  parameters  for  on-line  provision  of  maximum  manoeuvre  diagram 
data. 


It  is  still  in  its  original  condition  as  described  in  Reference  1. 

8.2.4  PS CON  Subroutine 


Subroutine  PSCON  is  called  from  the  main  program  when  10PT=1  (Option  4A) 
PSCON  reads  data  files  produced  bv  program  ACRAFT,  expressed  on  a  load  factor 
(or  turn  rate),  Mach  number  grid  for  each  pressure  height  requested,  and 
prepares  a  disk  file  P2«Al .CON  where  the  energy  parameter  is  expressed  on  a 
(Mach  number,  altitude!  grid  for  each  cf  the  requested  load  factor  or  turn  rate 
grid  prints.  P24A1 .CON  will  then  be  used  as  the  input  file  when  running  P4V2  t 
produce  contour  plots. 

In  addition  this  subroutine  now  calculates  the  lift  limit  boundarv  by 
working  out  for  each  Mach  number  the  height  at  which  the  load  factor  (or 
turnrate)  corresponding  to  the  maximum,  lift  coefficient  is  equal  to  the 
requested  load  factor  (or  turn  rate). 

This  subroutine  has  been  largely  rewritten.  The  modifications  are 
labelled  K2  ,K4  in  the  program  listing. 

8.2.5  PSDIFF  Subrout  me 


Subroutine  PSDIFF  is  called  from  the  main  program  when  IOPT=4(Opt ion  4D) 
PSDIFF  performs  the  functions  of  PSCON  with  two  separate  files  produced  by 
ACRAFT.  Subroutine  PSDIFF  produces  three  data  files,  P24A1.CON,  P24A2.CON  and 
P24DIF.CON  (one  for  each  aircraft  and  the  differential  data)  any  of  which  can 
then  be  processed  by  program  P4V2 . 

This  subroutine  has  only  been  slightly  changed  and  is  described  in 
Reference  1.  The  modifications  are  labelled  K2,K4  in  the  program  listing. 

8.2.6  RATE!  Subroutine 

Subroutine  RATE1  is  called  from,  the  mam.  program  when  IOPT=  2  ( opt  ion  4B). 
RATE  1  reads  data  produced  by  program  ACRAFT  expressed  as  energy  parameter  verso 
turn  rate  (or  load  factor)  for  a  range  of  altitudes  (pressure  height  or  energy 
height)  and  Mach  numbers,  and  prepares  a  plotter  file  P2.PLT  ready  for  plotting 
A  Separate  plotter  page  is  produced  for  each  value  of  the  altitude  variable. 

This  subroutine  has  been  slightly  changed  to  allow  for  the  choice  cf  loa 
factor  as  the  x  axis.  The  modifications  are  labelled  K4  in  the  program  listing 
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8.2.7  SEPCON  Subroutine 


Subroutine  SEPCON  is  called  from  the  main  program  when  10PT=8 (opt  ion  4C-)  . 
SEPCON  reads  data  files  produced  by  program  ACRAFT,  expressed  on  a  turn  rate  (or 
load  factor),  Mach  number  grid  for  each  pressure  heights  requested  and  prepares 
a  file  P 240 .CON  where  the  energv  parameter  is  expressed  on  a  Mach  number,  turn 
rate  (or  load  factor)  grid  for  each  of  the  requested  heights.  P24C-.C0N  will  then 
he  input  to  program  P4V2  to  produce  contour  plots. 

In  addition  this  subroutine  also  transmits  to  P24G.CON  the  lift  limit 
boundaries  calculated  in  program  ACRAPT. 

This  subroutine  was  entirely  rewritten  for  Version  2. 

8.2.8  RATE 2  Subroutine 


Subroutine  RATF2  is  called  from  the  main  program  when  I  OPT=  3(  opt  l  or.  4C). 

RATE2  reads  data  produced  by  program  ACRAPT  representing  the  optimum 
energv  parameter  and  correspond ing  turr.  rate  (or  load  factor)  for  a  range  of 
energv  states  and  load  factors  (or  turn  rates’1.  RATE2  then  prepares  the  plotter 
file  P2.PLT  ready  for  plotting. 

This  subroutine  has  been  entirely  rewritten  for  Version  2.  Previously 
(Version  1)  it  was  necessary  to  run  P4  to  obtain  these  plots.  It  is  now 
possible  (and  more  accurate)  to  create  the  plotter  file  directly  from  P2V2 
because  of  the  alternative  turn  rate  grid  variable. 

8.2.9  MMP  Subroutine 


Subroutine  MMP  is  called  by  the  main  program  when  10PT=6  (Option  4F). 

MMP  is  a  conversational  routine  to  assist  in  calculating  data  for  the 
Maximum  Manoeuvre  Persistance  Diagram.  Fuel  available  for  optimum  manoeuvres  at 
discrete  range  intervals,  together  with  energy  state,  maximum  sustained  turn 
rate  and  fuel  flow,  are  requested  as  terminal  input,  and  the  number  of  sustained 
turns  at  each  range  is  provided  as  terminal  output. 

This  subroutine  has  not  been  modified  in  Version  2  and  is  described  fully 
lr.  Reference  1  . 

8.2.10  GPIP  Subroutine 


Subroutine  GRID  is  called  by  the  main  program  when  10PT  =  7  (Option  42). 

GRID  read*  the  data  files  produced  by  program  ACRAFT,  expressed  on  a  Mac*', 
number,  load  factor  (or  turn  rate)  grid  for  each  height  requested  and  prepares  a 
d l si  file  P2.PLT  for  plotting.  Turn  rate  (or  load  factor)  and  energy  parameter 
are  plotted  against  Mach  number  for  each  load  factor  (or  turn  rate)  requested. 

A  separate  page  of  plots  is  produced  for  each  value  of  the  height  variable 
(energy  height  or  pressure  height). 

These  plots  are  intended  as  a-  overview  of  the  data  grid  and  are  produced  with  a 
predetermined  scale. 

This  subroutine  has  he« -  modified  to  allow  for  the  additional  optic*  of 
seler’i'g  load  factor  or  tjr  ra'e  a"  one  of  the  grid  variables.  The 
label  led  V-  : *  the  program  listing. 


ns  are 


f 

8.3  Subroutine  Library  P24LIB 

PI^LIF  contains  a  number  of  small  utility  subroutines.  Most  of  these 
subroutines  have  been  slightlv  modified  in  Version  2  mainly  to  cater  for  the 
extra  information  being  presented.  The  mod  1 f i c a t  i  on s  are  labelled  K3  and  Kb  in 
the  program  listing. 

8.3.1  INTRP  and  BADINP  Subroutines 

Routine  INTRP  (numerical  interpolation)  and  routine  BAPINP  (input  error 
text'  are  obtained  from  the  file  P1V2  and  are  described  in  Chapter  7, 

8.3.2  I  NLA?  Subroutine 


INL.Af  is  called  from  subroutine  PSDIFF  to  read  an  alphanumeric 
identifying  header  from,  an  input  disk  file  produced  by  program  ACRAFT. 

8.3.3  PL7LAE  Subroutine 


PLTLAB  outputs  the  identifying  header  read  by  INLAB  to  a  plotter  file. 

8.3.4  UNITS  Subroutine 


UNITS  is  called  by  P2IN  to  allocate  text  strings  to  variables  for  use  in 
terminal  dialogue  and  plotter  axis  labelling,  dependent  on  the  unit  system  and 
the  energy  variable  selected. 

8.3.5  VRLAB  Subroutine 


WRLAE  writes  the  descriptive  header  read  by  INLAB  onto  an  output  disk 
file  to  be  used  as  input  program  P4V2. 


9.  PROGRAM  P4V2  DESCRIPTION 
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> 


9.1  Loading  c f  Frcgram  on  DEC  Svstem.-lO 

The  cvtwr  ands  for  loading  program  PAV2  and  saving  a  core  image  on  d.s'e  or. 
DEC  s  v  s  t  e  tn-  1  0  are: 

.  LOAD  PAY2  ,  PALI  F,.  SEARCH,  P24  LIB -'SEARCH 
SAVE  PAV2 

The  routine  contained  in  these  files  are  listed  in  the  following  table 
and  described  ir  this  chapter. 


p-v: 

MA I  N 

p_n?_A  F-05C.C 

PCON  PES 

01  TXT 

PLOTD  GRIDC 

BOUND  ETSEND 

P-LIB 

COST 

DIAG  SMOOTH 

PROMPT  REALIN  DASH 

F2  a  LI  B 

INTRP 

BADINP  INLAB 

PLTI.AB  UNITS 

9.2  File  PAV2 

Most  routines  contained  in  this  file  have  been  restructured  and 
rewritten.  Their  purpose  is  described  here  and  further  detail1-  are  available  in 
the  source  program  listing  which  is  extensively  documented. 

Q  .  2  .  1  PA  m_A1N  Program 

The  program  PAY2  processes  data  files  produced  by  program  F2V2  and 
produces  contour  plots  for  off  line  plotting.  The  contour  plots  or  file  F-.PLT 
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9.2.2  PAP?- A  Subroutine 


Subroutine  P-PS-A  is  called  from  the  main  program,  when  options  AA  or  AD 
are  selected.  I:  us*-  t  h>-  input  data  filename  to  differentiate  between  requests 
for  simple  air  raft  or  comparison  plot*. 


1 


P4PS4A  controls  the  conversational  dialogue  which  allows  the  user  to  vary 
several  features  of  the  plots:  scale  of  plots,  number  and  level  of  contours, 
selection  of  load  factors  (or  turn  rate),  smoothness  of  plots,  inclusion  of 
energv  state  contours. 

P4PS4A  also  controls  the  setting  up  and  plotting  of  identifying  text  and 
axes,  the  reading  of  grid  data  and  boundary  points  from  file  P24A.CON,  the 
plotting  of  energy  rate  contours  together  with  the  lift  limit  boundary  and  the 
plotting  of  the  superimposed  energy  state  contours 

9.2.3  P4PS4G  Subroutine 


Subroutine  P4PS4G  is  called  from  the  main  program  when  option  4G  is 

chosen . 


P^PS4G  controls  the  conversational  dialogue  which  allows  the  user  to  vary 
several  features  of  the  plots:  scale  of  plots,  number  and  level  of  contours, 
selection  of  heights  (or  energy  states),  smoothness  of  plots,  inclusion  of  load 
factor  and  turn  radius  grid  (if  the  grid  selected  i6  turn  rate). 

P4PS4G  also  controls  the  reading  of  grid  data  and  boundary  points  from 
file  P24G.CON,  the  plotting  of  energy  rate  contours  together  with  the  lift  limit 
boundary,  and  the  plotting  of  the  superimposed  constant  load  factor  lines  and 
structural  load  factor  boundary  when  turn  rate  has  been  selected  as  a  grid 
variable.  Only  the  points  within  the  lift  limit  boundary  are  plotted  whereas 
contours  outside  the  structural  limit  boundaries  are  drawn  with  a  dashed  line. 

9.2.4  BOUND  Subroutine 


Subroutine  BOUND  is  called  by  subroutine  P4PS4A  or  P4PS4G  to  plot  the 
lift  limit  boundary  line  and  label  the  end  of  the  line. 

9.2.5  ETSBND  Subroutine 


Subroutine  ETSBND  is  called  by  subroutine  P4PS4G  to  plot  and  label  the 
line  corresponding  to  the  structural  load  factor  limit  on  a  (Mach  number,  turn 
rate)  grid. 

9.2.6  GRIDC  Subroutine 


Subroutine  GRIDC  is  called  by  subroutine  P4PS4G.  It  controls  the 
plotting  of  lines  of  constant  load  factor  and  lines  of  constant  turn  radius  on  a 
(Mach  number,  turn  rate)  grid. 

9.2.7  PCON  Subroutine 

Subroutine  PCON  is  called  by  subroutines  P4PS4A  or  P4PS4G  to  control 
calls  to  the  software  routine  PLOT.  Plotting  is  achieved  by  calls  to  PLOTD 
(P4V2)  for  straight  lines  between  points  or  calls  to  SMOOTH  (P4LIB)  for 
cubic  arcs  between  points. 

9.2.8  PES  Subroutine 

Subroutine  PER  is  called  by  subroutine  P4PS4A  to  control  the  plotting  of 
Es  contours  by  calling  either  PLOTD  or  SMOOTH  for  straight  line  or  smoothed 
plots,  respectively. 
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9.2.9  PLOT!)  Subroutine 


PLOTD  is  a  interface  routine  for  the  CALCOMP  software  routine  PLOT.  This 
routine  has  been  transferee!  from  Version  1. 

9.2.10  Ol'TXT  Subroutine 


Subroutine  Ol'TXT  is  called  by  PCON  at  the  start  of  a  contour  line  to 
output  the  contour  level  on  the  plotter. 

This  subroutine  was  derived  from  Version  1. 

*5.3  Subroutine  Library  P4L1R 


Most  subroutines  put  together  in  P4LIB  were  taken  from  a  set  of  general 
purpose  library  routines  available  to  all  ARL  staff.  These  routines  are 
outlined  bneflv  below  but  mure  details  are  available  in  Reference  1. 

9.3.1  CONT  Subroutine 


Subroutine  CONT  is  called  by  P4PS4A  or  P4PS4G  to  draw  contour  lines  for  a 
variable  stored  in  a  two  dimensional  array  (i.e.  PS(M,H)). 

9.3.2  DIAG  Subroutine 


DIAG  is  a  routine  called  by  CONT  to  interpolate  along  a  mesh  diagonal 
using  a  simple  biquadratic  representation  of  the  mesh. 

9.3.3  SMOOTH  Subroutine 


SMOOTH  is  a  routine  which  smooths  the  contour  plots  by  joining  successive 
points  on  the  contour  with  cubic  arcs  having  tangential  coincidence  at  their 
common  points  (knots).  The  tangential  slope  is  set  equal  to  the  mean  of  the 
slopes  of  the  straight  lines  that  would  otherwise  join  points  on  either  side  of 
the  knot.  When  smoothing  is  selected,  an  additional  call  to  SMOOTH  is  required 
to  complete  the  current  contour  before  proceeding  to  the  next  contour  request. 

9.3.4  PROMPT  Subroutine 

PROMPT  is  a  routine  called  by  SMOOTH  and  REALIN  which  writes  a  text  cm 
the  user's  terminal. 

9.3.5  REALIN  Subroutine 


REALIN  is  a  routine  called  by  SMOOTH  to  accept  real  numerical  data  from 
the  user's  terminal. 

9.3.6  DASH  Subrout  me 

DASH  is  a  new  subroutine  called  by  subroutine  PLOTD  (P4V2)  to  plot  a 
dashed  line  between  two  points. 

0 .4  Subroutine  Library  P24LIB 

This  library  has  been  described  in  Paragraph  8.3. 
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10,  EXAMPLES 


The  previous  chapters  have  discussed  the  three  programs  in  the  suite. 

This  chapter  presents  examples  of  the  interactive  dialogue  used  when  running  the 
programs,  and  of  the  tabular  outputs  and  plots  produced  when  using  the  programs 
in  a  co-ordinated  fashion.  Any  given  aircraft  is  represented  by  files  of 
propulsion,  aerodynamic  and  configuration  data,  together  with  the  source  and 
core  image  of  the  FORTRAN  programs.  A  suitable  naming  convention  is  required  to 
provide  efficient  housekeeping  of  these  files.  These  files  can  be  stored  on  DEC 
tape,  since  changes  to  one  or  more  files  will  be  needed  as  data  banks  and 
programs  are  developed. 

10.1  Opt i or  4 A 

This  option  corresponds  to  the  options  1,6,  15  and  20  in  the  overview  on 
Figure  2.2. 

Figure  10.1  shows  a  sample  dialogue  when  running  the  three  programs 
ACRAFT,  P2V2,  P4V2  to  produce  a  set  of  SEP  (Ps)  contours  plots  for  specified 
load  factors  (option  1  in  Figure  2.2).  Program  ACRAFT  produces  a  file 
(EXAMPL.TXT)  for  printing  and  a  file  (EXAMPL  .Nl’M)  for  input  to  program  P2V2. 

Figure  10.2  shows  a  portion  of  the  file  EXAMPL.TXT  and  illustrates  the 
main  features  of  the  listed  output.  The  header  page  identifies  the  run  and 
echoes  the  configuration  data  supplied  in  the  configuration  file.  In  the 
subsequent  pages  energy  rate  and  turn  rate  data  are  tabulated  at  each  grid 
point.  Also  listed  for  each  Mach  number  are  the  turn  rate  and  load  factor 
corresponding  to  the  lift  limit  boundary.  Energy  rate  values  outside  the  lift 
limited  envelope  are  flagged  with  the  values  of  -9999.99. 

The  Ig  load  factor  SEP  contour  plot  produced  using  the  dialogue  in  Figure 

10.1  is  shown  in  Figure  10.3.  Contours  are  labelled  with  a  number  indicating 
the  value  of  the  parameter  plotted,  as  shown  on  the  right  hand  side  of  the  plot. 
The  zero  SEP  contour  is  emphasised  by  a  dashed  line.  Superimposed  on  the  plot 
are  contours  of  constant  energy  state. 

Figure  10.4  illustrates  a  similar  dialogue  to  obtain  contours  of 
P s V / v f  for  specified  turn  rates  (option  6  in  Figure  2.2)  with  an  extract  of 
the  corresponding  output  shown  in  Figure  10.5  and  a  typical  plot  shown  on 
Figure  10.6.  4 

10.2  Opt  ion  4B 


This  option  corresponds  to  the  options  2,  7,  16,  21  in  the  overview  of 
Fi g  ire  2,2. 

Figure  10.7  shows  a  typical  dialogue  for  producing  plots  of  SEP  versus 
turn  rate  for  lines  of  constant  Mach  number  using  the  two  programs  ACRAFT  and 
P2V2  (option  2  in  Figure  2.2).  A  load  factor  increment  of  0.5  was  chosen  to 
provide  a  smooth  curve  when  the  values  are  plotted.  The  Mach  number  and 
altitude  selected  for  the  plots  are  nominated  when  running  P2V2 . 

Figure  10.8  shows  an  example  of  such  a  plot  where  each  curve  plotted  is 
identified  with  Mach  number  and  energv  state  annotation. 


149] 


10.3  Option  4C 

This  option  corresponds  to  the  options  11  and  13  in  the  overview  of 
Figure  2.2. 

Figure  10. o  shows  a  typical  dialogue  when  running  the  two  programs  ACRAFT 
and  P2V2  for  producing  a  maximum  manoeuvre  diagram  where  Ps  is  optimised  with 
respect  to  Mach  number  and  plotted  against  energy  state  Es  for  selected  values 
of  load  factor. 


An  extract  of  the  output  file  produced  is  shown  in  Figure  10.10  and  the 
corresponding  plot  in  Figure  10.11.  To  understand  the  meaning  of  the  data 
produced  a  summary  of  the  theory  used  may  be  useful. 

By  definition  r  _  G~s 
'  -  ~  dt 

therefore  the  time  to  climb  from  one  energy  state  Esj  to  another  ec -  if: 


■? 


d E 


•'Esl 

If  we  look  at  a  lg  Pg  contour  plots  (Figure  10.3)  it  can  be  shown 
theoret l ca 1 1 y  (Reference  4)  that  the  optimum  climb  for  minimum  time  will  be  the 
locus  of  the  points  at  which  the  lines  of  constant  Es  and  the  lines  of 
constant  Ps  have  common  tangents.  On  the  output  file  shown  in  Figure  10.10 
the  lg  load  factor  points  provide  the  profile  for  a  minimum  time  climb  which 
could  be  drawn  on  Figure  10.3.  The  output  file  also  provides  for  each  energy 
state  and  load  factor  the  fuel  flow  corresponding  to  the  Mach  number  and 
altitude  which  optimises  Ps.  The  right  hand  side  load  factor  column  provides 
optimum  sustained  turn  conditions  which  are  used  in  maximum  manoeuvre 
persistence  calculations.  Also  indicated  at  each  energy  state  are  the  time  to 
climb,  fuel  used  in  the  climb  and  range  obtained  by  integration  from  the  initial 
energy  state,  (More  details  are  available  in  Reference  1). 


Figures  10.12,  10.13  and  10.14  show  a  similar  run  but  with  lines  of 
constant  turn  rate.  In  a  similar  fashion  the  amount  of  fuel  used  to  climb  from 
energy  state  Eq.  to  another  E=.,  is: 


—  dE, 


and  option  4 C  can  be  run  with  Ps/wf  as  the  energy  parameter  to  provide  the 
optimum  climb  for  minimum  fuel  used,  which  can  then  be  used  in  conjunction  with 
lg  contours  of  Fs/wf  plotted  using  option  4A. 


Similarly  option  4C  can  be  run  with  the  two  other  energy  parameters 
PgV'wj  and  Ps/V  to  . rovide  optimum  climb  for  maximum  range  per  unit  mass 
of  fuel  and  maximum  range  respectively. 


10.4  Option  4D 

This  option  corresponds  to  the  options  3,  8,  17  and  22  in  the  overview  of 
Figure  2.2. 
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Figure  10.15  shows  part  of  the  dialogue  used  when  running  P2V2  and  P4V2 
for  producing  lg  SEP  differential  contours.  First  the  program  ACRAFT  written 
for  the  first  aircraft  was  run  using  data  from  the  store  file  AC1.STR  generating 
AC1.NUM  for  input  to  P2V2.  In  a  similar  manner  AC2.NUM  was  created  for  the 
second  aircraft.  This  part  of  the  dialogue  was  omitted  as  it  is  the  same  as  for 
option  4A . 

Figure  10.16  shows  the  plot  obtained  using  the  dialogue  above.  On  this 
plot  the  lift  limit  boundaries  for  both  aircraft  have  been  plotted  by  P4V2 
together  with  the  differential  Ps  contours. 

10.5  Option  4G 


This  option  corresponds  to  the  options  4,  9,  18,  23  in  the  overview  of 
Figure  2.2. 

Figure  1 h .  1 7  shows  a  typical  dialogue  used  when  running  the  three 
programs  ACRAFT,  t'2V2  and  P4V2  to  produce  a  set  of  SEP  contour  plots  at  a 
constant  altitude  on  a  Mach  number,  turn  rate  grid. 

Program  ACRAFT  produces  a  file  EXAMPL.NUM  for  input  to  program  P4V2  and  a 
file  EXAMPL.TXT  ready  for  printing  and  shown  in  Figure  10.18  which  is  similar  to 
the  one  described  in  Paragraph  10.2. 

A  plot  produced  using  the  dialogue  above  is  shown  in  Figure  10.19. 
Contours  are  labelled  at  their  start  with  a  level  number  which  is  identified  at 
the  right  hand  side  of  the  plot.  Superimposed  are  the  lines  of  constant  load 
factor  and  turn  radius.  The  lift  limit  and  structural  load  factor  boundaries 
are  also  generated  by  P4V2  with  the  contours  outside  the  structural  limit  drawn 
with  a  dashed  line  and  the  points  outside  the  lift  limit  boundary  ignored. 

Figure  10.20  shows  SEP  contour  plots  at  a  constant  altitude  but  on  a  Mach 
number,  load  factor  grid;  the  presentation  is  slightly  different  but  the  lift 
limit  has  again  been  generated  by  P4V2. 

10.6  Option  4Z 


This  option  corresponds  to  the  options  5,  10,  19,  24  in  the  overview  of 
Figure  2.2. 


Fig  10.21  shows  a  typical  dialogue  for  producing  a  plot  of  SEP  versus 
Mach  number  for  several  load  factors  at  a  given  altitude  using  the  two  programs 
ACRAFT  and  P 2 V 2 ;  the  plot  generated  is  shown  in  Figure  10.22. 


Plots  on  this  format  were  previously  used  (Reference  1)  to  obtain 
contours  of  Ps  on  a  Mach  number,  turn  rate  grid  at  a  specified  height;  this 
is  now  done  using  option  4G  and  this  format  does  not  present  much  interest  any¬ 
more  . 
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.RUN  ACRAFT 


FIG. 


MANEUVERABILITY  GRID  CALCULATION 

I  PIPER  I  PL  (0)  OR  i.  I.  (I)  UNITS  o 

AIRCRAFT  DATA  FILENAME  *"*  EXPMPL.STR 

PIRCRPFT  EXAMPLE  DATE  21-Jun-03  TIME  13:56  33.7 
3  AIM-9  -  GUNS  -  50-  FUEL 

PE.  PS/ UF.  P'hmxj/ UF  OR  P9/V  (I.  2.  3  OR  4)  7  1 

IS  HP  THE  HEIGHT  VARIABLE  -  Y 

PRESET  GRID  (Y.  N  OR  C/R)  7  n 

HPO  ,  HP  STEP  (FT),  NO-  OF  STEPS  7  o. 5000. 13 

MACHO*  MPCH  STEP  .  NQ.  OF  STEPS  '**  O.  .05.21 

GRID  FOR  LOPD  FPCTCR  (01.  OR  TURNRA7E  (1)  TO 

Go  .  G  STEP  .  NO.  QF  STEPS  '•  1,2.3 


FQUER  (MIL-1  00 .  MP  X 

-200  > 

■'  1 00 

UING  SUE 

EP  (IF  UOR 

TABLE) 

o 

NOZZLE  DEFLECTION 

( : F  UPR I PLlE ) 

P^MO'r  C'HERE  *  DEV  I A 

TICn 

-  ICPC.O 

OUTPUT  ( 

1  -.TEXT  ,  2-NQS.  .  3- 

BOTH) 

T*  •» 

3/P  FILENAME  FOP  UNIT  * 

■■  EXAMPL.TXT 

C/c‘  FILE 

NQME  FDR  UNIT  0 

"•  EXPMPL.  NUM 

cslculpt 

ION  : 

ALTITUDE 

HP  - 

0 . 0 

PT 

PLTIYJDE 

HP  - 

SOuO.  o 

F~~ 

altitude 

Up  a 

luouo .  u 

FT 

ALTITUDE 

HP  - 

15000. O 

FT 

altitude 

HP  - 

20AOO.  i.) 

FT 

ALTITUDE 

HP  - 

25<  n  >0.  t^i 

FT 

ALTITUDE 

HP  - 

20000  .  0 

PT 

ALTITUDE 

HP  - 

33(300 .  O 

PT 

altitude 

HP  - 

40000. 0 

P7 

altitude 

HP  - 

45O00.  r» 

FT 

ALTITUDE 

HP  - 

50000 . 0 

FT 

altitude 

HP  - 

55000.  O 

F~ 

ALT 1TUDE 

HP  - 

COOOU .  0 

FT 

STOP 

END  OF  £ 

XEC'JTICN 

CF-j  TIME 

»  13.05  EL 

PPSED 

TIME'-  2*4 

0 .  32' 

Zx:  T 

10.1(a)  Sample  Dialogue  for  Option  4A 

(Load  factor  grid,  Ps  as  energy  parameter) 
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Sample  Output  Listing  of  Program  ACRAt'T:  Header  Page 

(Load  factor  grid,  P  as  energy  parameter) 


ALTITUDE  HP  — 


0.0  FT 


FIG. 


MAX  MACH  «  1.0 


TABULATED 

VALUES 

:  PS 

FT/S 

TURNRATE 

- 

DEG/S 

I  MACH 

TAS 

ES 

H 

KF 

NO 

KT 

FT 

FT 

L3/HF 

2  .050 

33.1 

48. 

0. 

7418. 

3  .100 

66.1 

194. 

0. 

7479. 

A  .150 

99.2 

436. 

0. 

7541  . 

5  .200 

132.3 

775. 

0. 

7603. 

6  .250 

165.4 

1211. 

0. 

7664. 

7  .300 

198.4 

1743. 

0. 

7727. 

8  .350 

231.5 

2373. 

0. 

7795. 

9  .400 

264.6 

3099. 

0. 

7870. 

10  .450 

297.7 

3923. 

0. 

7953. 

11  .500 

330.7 

484  3. 

0. 

8044. 

12  .550 

363.3 

5860. 

0. 

3139. 

L»J 

O' 

o 

o 

396.9 

6973. 

0. 

8236. 

14  .650 

430.0 

3184. 

0. 

9333. 

15  .700 

463.0 

9*91. 

0. 

8422. 

16  .750 

496. 1 

10896. 

0. 

8496. 

17  .800 

529.2 

12397. 

0. 

85*9. 

18  .350 

562.3 

13995. 

0. 

3576. 

19  .900 

595.3 

15690. 

0. 

8580. 

20  .950 

629.4 

17482. 

0. 

3569. 

211.000 

661.5 

19370. 

0. 

8550. 

TRMAX  EMAX 
DEC/S  G  UNITS 

1.00  3.00  5.00 


0.00 

0.22 

-48.5 

-9999.9 

-9999.9 

0.00 

93.40 

161 .78 

0.00 

0.43 

1.8 

-87''.  3 

-3204.0 

0.00 

46.70 

80.89 

0.00 

0.77 

40.4 

-412.7 

-1486.8 

0.00 

31.13 

53.93 

6.28 

1.26 

80.2 

-286.4 

-1046.8 

0.00 

23.35 

40.45 

10.52 

1.38 

117.5 

-189.6 

-814.9 

0.00 

18.68 

32.36 

13.47 

2.64 

142. 1 

-107 . 1 

-642.1 

0.00 

15.57 

26.96 

16.05 

3.55 

162.7 

-33.5 

-500.7 

0.00 

13.34 

23.11 

18.48 

4.59 

179.1 

34.8 

-379.5 

0.00 

11.68 

20.22 

20.84 

5.77 

191.3 

83.7 

-271.4 

0.00 

10.38 

17.98 

23. 17 

7.09 

198.3 

101 .5 

-171 .7 

0.00 

9.34 

16,18 

25.48 

8.55 

200.9 

112.1 

-77 . 9 

0.00 

3.49 

14.71 

27.78 

10.  14 

196.' 

114. 1 

-4».  3 

0.00 

'.78 

13,  m? 

30.07 

11.88 

185.0 

100.2 

-42.6 

0.00 

7.18 

12.44 

32.36 

13.75 

160.0 

37.6 

-57.2 

0.00 

6 . 67 

11.56 

34.64 

15.77 

120.9 

50.8 

-89 .  ? 

0.00 

6.23 

10.79 

36.93 

17.92 

63.4 

-5.5 

-143.3 

0.00 

5.8“ 

10. 1 1 

39.21 

C\J 

O 

-38.4 

-109.9 

-252.9 

0.00 

5.“9 

9.52 

41.50 

22.6* 

-219.9 

-299.2 

-457.7 

0.00 

5.19 

8.09 

43.78 

25.21 

-5*6.8 

-641.1 

-829.3 

0.00 

4.92 

8.51 

46.07 

27.92 

-979  .7 

-1088.9 

-1707.2 

0.00 

4.67 

8.09 

LOAD  FACTR 


10.2(b)  Sample  Output  Listing  of  Program  ACRAFT:  Extract 
(Load  factor  grid,  Ps  as  energy  parameter) 
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RUN  ACRAF7 


MANEUVERABILITY  GRID  CALCULATION 

IMPERIAL  (0)  OR  S.  !  .  (1)  UNITS  ?  0 

AIRCRAFT  DATA  FILENAME  7  EXAMPL.STR 

AIRCRAFT  EXAMPLE  DATE  11 -Oct-63  TIME  1  5 « 27  06.4 
2  RlH-9  ♦  GUNS  -  SOX  FUEL 

PS.  PS/UF.  PS*V/UF  OR  PS/V  (1.  2,  3  OR  4 )  ?  3 

IS  HP  THE  HEIGHT  VARIABLE  7  Y 

PRESET  GRID  (Y,  N  OR  C/R)  7  N 

MPO  .  HP  STEP  (FT).  NO.  OF  STEPS  7  0.5000.13 

MACHO.  MACH  STEP  .  NO.  CF  STEPS  7  0.. 025. 41 

GRID  FOR  LOAD  FACTOR  (0).  OR  TURNRATE  (1)  7 « 

RO  .  TR  STEP  .  NO.  OF  STEP’S  7  0.4,3 

F’OUER  (MIL- 100.  MAX-200)  7  100 

WING  SWEEP  (IF  VARIABLE)  °  0 

NOZZLE  DEFLECTION  (IF  VAR I  ABLE  >  7  o 

ATMOSPHERE  .  DEVIATION  7  ICAO.O 

OUTPUT  (  1 -TEXT  .  2— NOS,  .3“ BOTH)  7  3 

3/P  FILENAME  FOR  UNIT  6  7  EXAMPL.TXT 

C/c’  CILENAME  FOR  UNIT  6  7  EXAMPL.NUM 

CALCULATION  « 


ALTITUDE 

HP 

- 

0.0 

"7 

ALTITUDE 

HP 

- 

5000. 0 

FT 

ALTITUDE 

HP 

- 

1  OC'OO .  0 

FT 

altitude 

HP 

- 

1 5000 . 0 

FT 

ALTITUDE 

HP 

- 

20000 . 0 

CT 

altitude 

HP 

- 

25000.0 

PT 

ALTITUDE 

HP 

- 

30000. 0 

er 

altitude 

HP 

- 

35000. 0 

FT 

altitude 

HP 

- 

40000 . 0 

FT 

altitude 

HP 

- 

45000. 0 

FT 

altitude 

HP 

- 

50000 . 0 

FT 

altitude 

HP 

- 

55000. 0 

FT 

alt  I TUDE 

HP 

- 

60000 . 0 

FT 

STOP 

END  OF  EXECUTION 

CPU  TIME'  40.12  ELAPSED  T  I  ME :  S *  37 . 32 
EXIT 


FIG.  10.4(a)  Sample  Dialogue  for  Option  4A 

(Turn  rate  grid,  PsV/wf  as  energy  parameter) 


. RUN  P2V2 


COMDAT  PERFORMANCE  PROCESSING 
OPTION  OR  (CR)  FOR  HELP  '  4A 

4  A  PS  CONTOURS  ON  ALT  US  MACH  GRID  (P24A1.C0N) 

DRTR  BRSE  FILENAME  ■  EXAMPL.NUM 

mmm  DRTR  RRE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PSU/ 1 OOOUF  *** 
.*»  CONSTANT  INC.  OF  TURN  RATE  *** 

REPLY  “YES".  "NO".  "ALL"  OR  “END"  i 
OUTPUT  ORTA  FOR  TR  -  0.00  7  ALL 

END  OF  EXECUTION 

CPU  TIME>  12.73  ELAPSED  TIME ‘  3>2.3& 

EXIT 

.RUN  P4UD 


INPUT  FILENAME  :  P24A1.C0N 

mmm  ORTA  ARE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PSU/1000WF  ««* 
mmm  CONSTANT  INC.  FOR  TURN  RATE  «»* 

CONTOUR  PLOTTING 


SCALES  IN  UNITS/IN  OF 

PLOT  -  MACH. 

ALT  i  .2. 10000 

SMOOTHED 

AND  TEXT 

UPEO 

CONTOURS 

N 

pLQT  ES' 

CONTOURS 

-*>  Y 

CONTOUR 

levels  - 

START 

.  STEP.  NO. 

«  5000.5000.15 

PLOT  TR 

-  0 . 00 

'  M 

PLOT  TR 

-  4  .  oo 

-  - 

CONTOUR 

levels  - 

STAR" 

.  STEP.  NC. 

i  -300% 50-13 

»LOT  TP 

-  5 .  GO 

"  N 

STOP 

END  OF  EXECUTION 

CPU  TIME'  11.32  ELAPSED  T ! ME •  2 <4.32 
EX  !  T 


FIG.  10.4(b) 


Sample  Dialogue 


for  Option  4A 


(cont) 


(58) 


ALTITUDE  HP  =  tOOOO.O  FT 


MAX  MACH  = 


TABULATED  VALUES  :  PS*V/WF  -  FT»*2/LB. S/1000 
LOAD  FACTR  -  0 


I  MACH 

TAS 

ES 

H 

WF 

TRMAX 

EMAX 

NO 

iCT 

FT 

FT 

LB/ HR 

DEG/ 3 

0  UNITS 

0.00 

4.00 

8.00 

A  .075 

47.9 

10101 . 

10000. 

5543. 

0.00 

0.21 

-3.7 

1.00 

-4.0 

1.02 

-4.3 

1.06 

5  .100 

63.3 

.0180. 

10000. 

5575. 

0.00 

0.30 

-3.0 

1.00 

-3.4 

1.03 

-4.5 

1 .  10 

6  .125 

79.3 

10282. 

10000. 

5609. 

0.00 

0.41 

-2.2 

1.00 

-2.3 

1.04 

-4.7 

1.16 

7  .150 

95.7 

10406. 

10000. 

5643. 

0.00 

0.54 

-0.3 

1.00 

-1.6 

1.06 

-4.7 

1.22 

3  .175 

in  .7 

10552. 

10000. 

5679. 

0.00 

0.69 

1  .00 

-0.3 

1.08 

-4 . 4 
1.29 

9  .200 

127.7 

10722. 

10000. 

5716. 

0.00 

0.37 

3.6 

1 .00 

1  1 

1.10 

-3.Q 

1 . 37 

10  .225 

i  u  3 . 6 

10913. 

10000. 

5755. 

2.98 

1.07 

6 . 5 
1.00 

4  .  1 
1.1? 

-3.2 

1.45 

11  .250 

159 . 6 

11127. 

10000. 

5795. 

5.70 

1.30 

10.0 

1.00 

6.9* 

1 .  16 

-2.3 

1 .54 

12  .275 

175.5 

11364. 

10000. 

583.6. 

'.39 

1.55 

13.9 

1.00 

10.2 

1 .  19 

-1.1 

1.63 

13  .300 

191 .5 

11623. 

10000. 

5879. 

3. '2 

•.37 

18.0 

1 . 00 

13.3 

1.22 

0.2 

1.72 

l a  . 325 

207.5 

11905. 

1000c. 

5923. 

9.38 

21. a 
1.00 

17.3 

1 . 26 

l .  9 
1.32 

15  .350 

223.4 

12210. 

10000. 

5969. 

10.92 

2.45 

24.0 

1 . 00 

21  .? 

1 .29 

1.92 

16  .375 

239.4 

12537. 

10000. 

6017. 

11.91 

2.70 

28.5 

1.00 

24.3 

1 . 3? 

5.9 

2.02 

17  .400 

255.3 

12886. 

10000. 

6066 . 

12.35 

;.  17 

72 . 0 
1.CC 

27.? 

1. 37 

a  7 

2.  12 

13  .425 

271.3 

'3258. 

10000. 

5117. 

13.75 

3.56 

35.6 

1.00 

.  3 
1.41 

11.1 

2.22 

19  .450 

287.2 

13653. 

10000. 

6170. 

14.64 

3.98 

30.0 

1.00 

33.2 

1 .  45 

14.2 

2.33 

20  .475 

303.2 

14070. 

10000. 

6224. 

15.51 

4.42 

42.4 

1.00 

35.9 

1  .  4Q 

16.3 

2.44 

21  .500 

315.2 

14510. 

10000. 

6280. 

16.3' 

4.89 

45.6 

1.00 

?8.4 

1.54 

16. 9 
2.54 

22  .525 

335.' 

14972. 

10000. 

6338  . 

17.21 

5.38 

48.6 

1.00 

40.7 

1.58 

17.2 

2.65 

23  .550 

351.1 

15«57. 

10000. 

6398. 

18.06 

5.89 

51.4 

1.00 

42.3 

1.63 

17.1 

2.76 

24  .575 

367.0 

15964. 

10000. 

6459. 

18.30 

6.43 

53.9 

1.00 

Uu .  6 
1.63 

16.7 

2.87 

25  .600 

383.0 

16494 . 

10000. 

6523. 

19.73 

6.99 

56.0 

1.00 

46.0 

1 .72 

15 . 9 
2.08 

26  .625 

399.0 

17046. 

10000. 

6588. 

20.56 

7.57 

57.9 

1.00 

47  .  1 
1.77 

14.3 

3.09 

27  .650 

414.9 

17621. 

10000. 

5054 . 

21.33 

3.18 

59.0 

47.4 

12.5 

FIG.  10.5  Sample  Output  Listing  of  Program  ACRAFT 

(Turn  rate  grid,  P  V/wf  as  energy  parameter) 
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RUN  ACRAFT 


MANEUVERABILITY  GRID  CALCULATION 

IMPERIAL  (0)  QR  S.I.  (1)  UNITS  7  0 

AIRCRAFT  DATA  FILENAME  7  DUMMY. STR 

AIRCRAFT  EXAMPLE  DATE  U-0ct-03  TIME  16*05  56-9 

AIR  SUPERIORITY  i  2  AIMS  *  502  FUEL 

PS.  PS/UF.  PS-U/UF  OR  PS /V  (1.  2.  3  OR? 4)  71 

IS  HP  THE  HEIGHT  VARIABLE  7  Y 

PRESET  GRID  (Y,  N  OR  C/R)  7  N 

HPO  .  HP  STEP  (FT).  NG.  OF  STEPS  7  0.10000.2 

MACHO.  MACH  STEP  .  NO.  OF  STEPS  7  .4.. 2. 6 

GRID  FOR  LOAD  FACTOR  <0>.  OR  TURNRATE  (1)  70 

GO  .  G  STEP  .  NO.  OF  STEPS  7  i.,5.15 

POWER  (MlL-100. MAX-200)  7  200 

WING  SWEEP  (IF  UARIABLE)  7  0 

NOZZLE  DEFLECTION  (IF  UARIABLE)  7  Q  * 

ATMOSPHERE  .  DEVIATION  7  ICAQ.O 

OUTPUT  ( 1-TEXT. 2-NCS. . 3- BOTH)  7  3 

Q/P  FILENAME  FOR  UNIT  6  7  EXAMPL.TXT 

0/P  FILENAME  FOR  UNIT  5  7  EXAMPL.NUM 

CALCULATION  • 

ALTITUDE  HP  -  0.0  FT 

ALTITUDE  HP  -  10000.0  FT 

STOP 

END  OF  EXECUTION 

CPU  TIME*  4.43  ELAFSED  TIME*  3*^. 56 
EXIT 

. RUN  P2V2 

COMBAT  PERFORMANCE  PROCESSING 
OPTION  OR  <CR)  FOR  HELP  *  4B 

4B  PS  US  OMEGA/ETA  FOR  GIVEN  HEIGHT.  (P2.PLT) 

DATA  BASE  FILENAME  •  EXAMPL.NUM 

*•*«■  DATA  ARE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PS 
**«  CONSTANT  INC.  OF  LOAD  FACTOR  mmm 
SCALES  IN  UN  ITS/ INCH  OF  plot 

TURN  RATE/ LOAD  FACTOR  .  ENERG Y  VARIABLE  •  4,200 
MAXIMUM  FOR  ENERGY  AXIS  »  400 


ENERGY 

AXIS  LENGTH  (IN)  *  7 

REFLY  - 
HEIGHT 

YES 

".  -NO".  ■'  ALL  ’*  OR 

0. 0  FT  7  N 

HEIGHT 

- 

10000.0  FT  7  Y 

MACH 

- 

.  400  7  N 

MACH 

- 

.600  7  N 

MACH 

- 

. 000  ?  Y 

MACH 

-  1 

.000  7  Y 

MACH 

-  1 

. 200  7  E 

END  OF  EXECUTION 

C°'J  TIME*  4.63  ELAPSED  TIME*  2*26.  DO 
EXJ  T 


FIG.  10.7  Sample  Dialogue  for  Option  4B 
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ii  -Mpfii 


RUN  OCRAFT 


I'lAuruVERHtlLn  Y  GRID  CALCULATION 

I  PIPER  I  Al  «0)  OR  3.1.  (1»  UNITS  ?  '» 

AIRCRAFT  DATA  FILENAME  7*  EXAMF'L.  3TR 

AIRCRAFT  EXAMPLE  DATE  lO-0ct-03  T 1 ME  10*34 
2  AIM-3  *  GUNS  *  sot  FUEL 

PS.  PS/Ur.  PS-U/UF  QR  P3/V  (1.  3.  3  OR  4  i  "»  1 

IS  HP  THE  HEIGHT  UARIAfcLc  T  N 

maximum  maneuver  -  y 

FISO  .  ES  STEP  (FT).  NO.  OF  STEPS  7  SOW.  50*  **J.  1 3 
GRID  FOR  LOAD  FACTOR  lO>.  OR  TIJRNRATE  (1)  'O 
GO  .  G  STEP  .  NO.  OF  STEPS  "  1 « 1.0 

POUER  ( mil- l OO. MAX- sou) 

UilNG  SViEEP  (IF  UARIArCe) 


NOZZLE  DEFLECTION  (IF  VARIABLE)  T  o 


ATMOSPHERE  .  DEO I AT  I  ON 
OUTPUT  ( 1- TEXT . 2-NOS. . 3-BQThj 
c/p  filename  FOR  UNIT  t- 
0/P  FILENAME  FOR  UNIT  j 


'•  I  COO.  U 

T  EXPORT. 7X7 
**  EXPORT.  nUM. 


CALCULATION 

ENERGY 

.“TATE 

EC 

ENERGY 

STATE 

EC 

- 

ENERGY 

STATE 

c; 

- 

ENERGY 

state 

ES 

- 

ENERGY 

STATE 

ES 

- 

r  NFRGY 

STATE 

E3 

- 

ENERGY 

state 

EL 

- 

ENERGY 

state 

ES 

- 

ENERGY 

•:.rn  te¬ 

FL 

- 

ENERGY 

ES 

- 

ENERGY 

state 

ES 

- 

ENERGY 

STATE 

ES 

- 

ENERGY 

STOP 

STATE 

ES 

- 

5*  »  »J.  1 1  FT 
i <  h 'Hu.  u  p; 
ISi  »«■»!),  IJ  FT 
2«U.iOu.O  FT 
iSrmO.o  P- 
Ji  m’m  io .  •  •  ft 
35*5* *0.0  FT 

41 II  HIM  .A  FT 
4 5000. O  FT 

FT 

SS4.H»0.0  FT 
EOMOO.O  FT 
G5OU0.H  FT 


END  cr  EXECUTTON 

CPU  Tine*  40.39  ELAPSED  TIME*  5*14.60 
EXIT 


.  RUN 


P2V2 


COMBAT  PERFORMANCE  FT,CCES3ING 
OPTION  OR  CR  )  FOR  HELP  •  4C 

•1C  MAXIMUM  MANEUVER  DIAGRAM  -  HMD.  >  PC 

DATA  BASE  FILENAME  '  EXPORT. NUn 

-**  DATA  ARE  IN  IMPERIAL  UNI"’-,  ENCRGV  PARAMETER 
•••  CONSTANT  INC.  CF  LOAD  "ACTOR  «  — 

S  CALE  IN  UNIT/ INCH  OF  PLOT 
E  3  .PS  *  lot  mo.  COO 

rs  MIN  VALUE.  P3  AXIS  LENGTH  (  IN-;.)  i  -7IUI.S 

F  F  c  L  y  -YES”.  -NO".  “A| _ •  OR  ’  End  "  • 

G  -  i. A* i  ?  nLi_ 

FnD  or  EXECUTION 

CPU  TIME*  3 ,?U  EL APF ED  TIME*  l*i0.lo 
EC' 


.  PLT  ) 

I  L  P’S 


G.  10.9 


Sample  Dialogue  for  Option  4 
(Load  factor  grid) 


•••  ENERGY  STATE  ES  s 


5000.0  FT  **• 


G 

1.00 

2.00 

PS 

FT/S 

199.5* 

163.60 

OMEGA 

DEG/S 

0.000 

5.629 

MACH 

0.508 

0.506 

TAS 

KT 

336.044 

336.060 

HP 

FT 

1 . 

0. 

WF 

LB/HR 

8059. 

3059. 

TIME 

S 

0.00 

FUEL 

LB 

0. 

RANGE 

N.M. 

0.000 

••*  ENERGY 

STATE  ES 

VI 

1.00 

2.00 

PS 

FT/S 

166.91 

159.16 

OMEGA 

DEG/S 

0.000 

9.659 

MACH 

0.608 

0.620 

TAS 

:<T 

398. 166 

906.518 

HP 

FT 

2982. 

2689. 

WF 

LB/ HR 

7779. 

7859. 

TIME 

S 

25.88 

FUEL 

LB 

57. 

RANGE 

N.M. 

2.591 

ENERGY 

STATE  ES 

G 

1.00 

2.00 

PS 

FT/S 

167.19 

130.52 

OMEGA 

DEG/S 

0.000 

9.507 

MACH 

0.633 

0.651 

TAS 

KT 

907.931 

919.731 

HP 

FT 

7651. 

7201 . 

WF 

L3/HR 

7030. 

7155. 

TIME 

S 

59.12 

FUEL 

LB 

115. 

RANGE 

N.M. 

5.659 

3.00 

9.00 

5.00 

9.11 

103.70 

19.85 

■156.25 

0.00 

9. 192 

12.587 

15.921 

12.963 

0.508 

0.508 

0.508 

0.508 

336.069 

336.069 

336.070 

336.069 

0. 

0. 

0. 

0. 

8059. 

8059. 

8059. 

3059. 

10000.0  FT  ••• 


3.00 

9.00 

5.00 

9.39 

101.55 

30.60 

-57.98 

0.00 

7.373 

9.79“ 

12.096 

10.577 

0.638 

0.697 

0.679 

0.661 

919.010 

931.908 

999 . 179 

936.171 

2223. 

1792. 

1266. 

1576. 

7966. 

8081  . 

8189 . 

3113. 

15000.0  FT  *** 


3.00 

9.00 

5.00 

3.90 

71.95 

-3.17 

-106.97 

o.on 

7.097 

9.963 

11.578 

9.220 

0.67“ 

0.691 

0.712 

0.639 

“35.265 

997.023 

962.126 

995.817 

6613. 

6153. 

5596. 

6201 . 

73'“. 

7939. 

7589. 

7922. 

FIG.  10.10 


Extract  of  Sample  MMD  with  a  Load  Factor  Grid 
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E0IXJ 


[€  5] 

1 


RUN  ACRAFT 


r 


MANEUVERABILITY  GRID  CALCULATION 


IMPERIAL  (O)  OR  S.l,  (1)  UNITS  ?  O 


AIRCRAFT  DATA  FILENAME 


?  EXAMPL. STR 


AIRCRAFT  EXAMPLE  DATE  li-0ct-03  TIME  15*12  21.0 
2  AIfl-9  ♦  GUNS  ♦  SOX  FUEL 


PS*  PS/UF.  PS*V/UF  QR  PS/V  11.  2.  3  OR  4)  ?  L 


IS  HP  THE  HEIGHT  VARIABLE 


MAXIMUM  MANEUVER 


ESO  .  ES  STEP  (FT),  NO.  OF  STEPS  7  0.10000.6 
GRID  FOR  LOAD  FACTOR  (O) .  OR  TURNRATE  (1)  ?1 


RO  ,TR  STEP 


NO.  OF  STEPS  7  0.4.4 


POWER  miL-100. MAX-200) 


UING  SWEEP  (IF  VARIABLES 


NOZZLE  DEFLECTION  (IF  VARIABLE)  ?  0 
atmosphere  .  DEVIATION  ?  IC 


OUTPUT  C 1-TEXT. 2-NOS. . S-BOTH) 


O/P  FILENAME  FOR  UNIT  6 


0/p  FILENAME  for  UNIT  3 


CALCULATION  « 

ENERGY  STATE  ES  - 
ENERGY  STATE  ES  - 
ENERGY  STATE  ES  - 
ENERGY  STATE  ES  - 
ENERGY  STATE  E3  - 
ENERGY  STATE  ES  - 
ENERGY  STATE  ES  - 
ENERGY  STATE  ES  - 
STOP 


7  EXAMPL.TXT 


EXAMPL. NUM 


0.0  FT 
10000.0  FT 
20000 .  O  FT 
30000.0  FT 
40000.0  FT 
50000.0  FT 
60000.0  FT 
70000.0  FT 


END  OF  EXECUTION 

CPU  TIME •  25.67  ELAPSED  TIME*  3*3.72 
EXIT 

.RUN  P2V2 


CQMOAT  PERFORMANCE  PROCESSING 
OPTION  OR  (CR)  FOR  HELP  •  4C 

4 C  MAXIMUM  MANEUVER  0 1 AGRAM  -  MMD.  (P2.PLT) 

DATA  8ASE  FILENAME  *  EXAMPL. NUM 

DATA  ARE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PS 
***  CONSTANT  INC.  OF  TURN  RATE 
SCALE  IN  UN  IT/ INCH  OF  PLOT 
ES  .PS  •  10000.200 

PS  MIN  VALUE.  PS  AXIS  LENGTH (  INS)  *  -000.6 

REPLY  -YES**.  “N0“.  ••ALL"  OR  -  END M  • 

TR  -  0.00  7  ALL 

END  OF  EXECUTION 

CPU  TIME*  5.05  ELAPSED  TIME*  1*0.66 
EXIT 


FIG.  10.12  Sample  Dialogue  for  Option  4C 
(Turn  rate  grid) 
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1 


ENERGY  STATE  ES  *  10000.0  FT 


OMEGA 

DEG/S 

0.000 

9.000 

8.000 

12.000 

11.968 

PS 

FT/S 

186.91 

163.  99 

99.99 

-19.57 

0.00 

G 

1.00 

1.73 

2.87 

9.01 

3.86 

MACH 

0.608 

0.591 

0.563 

0.593 

0.595 

TAS 

KT 

398. 166 

386.357 

367.269 

353.352 

355.209 

HP 

FT 

2982. 

3392. 

9029. 

99  73. 

99  13. 

WF 

L3/HR 

7779. 

7675. 

7507. 

7386. 

7902. 

TIME 

S 

0.00 

FUEL 

LB 

0. 

RANGE 

N.M. 

0.000 

•••  ENERGY 

STATE  ES  = 

20000.0 

FT 

OMEGA 

DEG/S 

0.000 

9.000 

8.000 

12.000 

8.315 

PS 

FT/S 

199.88 

111.62 

in. 314 

-167.73 

0.00 

G 

1.00 

1.81 

3.05 

4.91 

3.20 

MACH 

0.666 

0.659 

0.624 

0.689 

0.629 

TAS 

KT 

921.725 

913.925 

393.545 

4?7 .415 

366.99° 

HP 

FT 

12126. 

12933. 

1R141; . 

1  l^3C . 

1°0 16 . 

WF 

L3/HR 

6316. 

6225. 

6011. 

64Q1  . 

609° . 

TIME 

S 

60.28 

fuel 

LB 

ne. 

RANGE 

N.M. 

6.693 

•••  ENERGY 

STATE  ES  = 

30000.0 

FT  111 

OMEGA 

DEG/S 

0.000 

9.000 

3.000 

12.000 

5.356 

?S 

FT/S 

97.78 

93.50 

-9“ .59 

-900.06 

0.00 

'J 

1.00 

1.37 

3.26 

5.59 

2.34 

MACH 

0.798 

0.705 

0.6°6 

0.305 

CM 
0 
r  - 

0 

TAS 

KT 

958.125 

330. 260 

929 . 109 

996.000 

923.203 

HP 

FT 

20709. 

21805. 

2203«. 

10109. 

21882. 

Hr 

LB/ HR 

5080. 

9767. 

9702. 

5591. 

3795. 

TIME 

S 

192.70 

FUEL 

LB 

298. 

RANGE 

N.M. 

16.659 

FIG.  10.13  Extract  of  Sample  MMD  with  Turn  Rate  Grid 
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(0IX)  S/ld  -  Sd  3iba  JLDd3N3 


FIG.  10.14  Sample  Maximum  Manoeuvre  Diagram  for  Specified 
Turn  Rate 


!  C6  ] 


.  RUN  PCU: 


COMBAT  PERFORMANCE  PROCESSING 
OPTION  OR  (CR)  FOR  HELP  •  40 

4D  PS  DIFFERENT.  PLOT ( P24 A i . CON. P24A2 . CON . P24D I F . CON > 
DATA  BASE  FILENAME  *  AC  1  .  NUP1 

*<**  DATA  ARE  IN  IMPERIAL  UNITS*  ENERGY  PARAMETER  IS  PS 
***  CONSTANT  INC.  OF  load  FACTOR  **w 
REPLY  “YES".  "NO",  "ALL"  OR  "END"  • 

OUTPUT  DATA  FOR  N  -  1 . 00  ?  Y 

COMPARISON  FILENAME  •  AC2.NUM 

END  OF  EXECUTION 

CPU  TIME*  13.19  ELAPSED  TIME*  1*47.32 
EXIT 

. RUN  P4U2 


input  filename  *  P24Dif.ccn 

**■*  DATA  ARE  IN  IMPERIAL  units*  ENERGY  PARAMETER  IS  PS 
«*W*  CONSTANT  InC.  FOR  LOAD  "ACTOR  mmm 

contour  plotting 

SCALES  in  UNITS/IN  OF  PLOT  -  MACH.  ALT  :  .2*10000 

SMOOTHED  AND  TEXTURED  CONTOURS  7  N 
PLOT  ES  CONTOURS  7  n 
PLOT  G  -  I . 00  7  y 

CONTOUR  LEUCLS  -  START.  STEP.  NO.  •  -100.20.9 
3  "OP 

END  OF  EXECUTION 

cpu  T  i  me  *  c.  ‘39  elapsed  time*  uss.eo 

EXIT 


FIG.  10.15  Sample  Dialogue  for  Option  4D 
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.RUN  ACRAFT 


(MANEUVERABILITY  GRID  CALCULATION 

IMPERIAL  (O)  OR  3.1.  (1)  UNITS  ?  O 

AIRCRAFT  DATA  FILENAME  ?  EXAMPL.  3TR 

AIRCRAFT  EXAMPLE  DATE  12-0ct-tJ3  TIME  14*S2  09.1 
2  AIM-9  ♦  GUNS  -  50X  FUEL 

PS.  PS/UF,  PS*V/UF  OR  PS/U  Cl.  2.  3  OR  4 )  P  l 

IS  HP  THE  HEIGHT  UAR I  ABLE  ?  Y 

PRESET  GRID  (Y,  N  OR  C/R)  N 

HPO  .  HP  STEP  (FT).  NO.  OF  STEPS  P  0.10000.2 

MACHO.  MACH  STEP  .  NO.  QF  STEPS  "•  0.  .025.41 

GRID  FOR  LOAD  FACTOR  (0).  OR  TURNRATE  (1) 

RO  . TR  STEP  .  NO.  OF  STEPS  -  O. .5.41 

POWER  (MIL- 100. MAX-200)  P  100 

UING  SLEEP  (IF  VARIABLE)  "*  0 

NOZZLE  DEFLECTION  (IF  UARIAfcLE?  °  O 

ATMOSPHERE  .  DEVIATION  -  ICPO.C 

OUTPUT  (  1 -TEXT.  2-NOS.  .  3-BC^H;  3 

0/  °  “ILENAME  -CP  UNIT  S  -  EXAMPL .  TXT 

O/P  FILENAME  FOP  UNIT  6  EXAMPL. NUM 

CALCULATION  * 

ALTITUDE  HP  -  0.0  - T 

ALTITUDE  HP  -  10000.0  f- 

STOP 

END  OF  EXECUTION 

CPU  TIME:  25.49  ELAPSED  TIME'  3«4?.lu 
EXIT 


FIG.  10.17(a)  Sample  Dialogue  for  Option  4 G 

(Turn  rate  grid,  Pg  as  energy  parameter) 
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. RUN  P2V2 


COMBAT  PERFORMANCE  PROCESSING 
OPTION  OR  (CR)  PGR  HELP  »  4G 

4G  PS  CONTOURS-  ON  OPEGA/LF  VS  MACH  GRID  (P24G.CON) 
DATA  EASE  FILENAME  »  EXAMPl.nuM 

*****  DATA  ARE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PS 
***  CONSTANT  INC.  OP  TURN  ROTE  *** 

REPLY  -YES"  «  "NO"  .  •‘ALL"  OR  "END’  t 
OUTPUT  DATA  POP  H  -  O.^ALL 

END  OF  EXECUTION 

CPU  TIME*  12.17  ELAPSED  TIME*  1*0.24 
EXIT 


.RUN  D4V2 


INPUT  FILENAME  *  P24G.CCN 

m»*  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS 
mmm  CONSTANT  INC.  FOR  TURN  RATE  *** 

CONTOUR  PLOTTING 

scales  in  unit: /in  of  f-lo^  -  mach.*r,l-  •  .2.4 

SMCCTUED  AND  TEXTURED  CONTOURS  n 

F'UC"  _OAD  "AC .  "URN  RAD.  GRID  “■  r 

6LOT  h  -  0.  ^  Y 

CONTOUR  LEVELS  -  START.  STEP,  nC.  •  -200. SO. 9 
PLO"  h  -  10000.  ^  n 

S'.  TCP 

END  CP  EXECUTION 

CPU  TIME:  ;7.3t  ELAPSED  TIME*  1*43.  '.>4 

ex:- 


FIG.  10.17(b) 


Sample  Dialogue  for  Option  4G 


(cant) 
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HLT1TUDK  ltl> 


FIG.  10.18  Sample  Output  Listing  of  Program  ACRAFT:  Extract 
(Turn  rate  grid,  P  as  energy  [jaraniuter) 


.RUN  ACRAFT 


MANEUVERABILITY  GRID  CALCULATION 

IMPERIAL  (0)  OR  S.I.  CU  UNITS  ?  O 

AIRCRAFT  DATA  FILENAME  7  DUMMY . STR 

AIRCRAFT  EXAMPLE  DATE  1 1 -Qct-03  TIME  16*0u  06.0 

AIR  SUPERIORITY  •  2  A I MO  *  SOX  FUEL 

PS.  PS/UF,  PS*u/uF  OR  PS/V  (I.  2.  3  OR  4 )  ° 

IS  HP  THE  HEIGHT  VARIABLE  7  Y 

PRESET  GRID  (Y.  N  OR  C/R)  •’  N 

HPO  .  HP  STEP  (FT),  NO.  OF  STEPS  7  U.lUUUU.2 

MACHO.  MACH  STEP  .  NO.  OF  STEPS  7  0.-05.40 

GRID  FOR  LOAD  FACTOR  (0).  OR  TURNRA7E  (1)  70 

Go  .  G  STEP  .  NO.  Of  STEFS  '*»  1.2.4 

POWER  (MIL- 100.  MAX-200)  7  200 

WING  SWEEP  (IF  VARIABLE)  7  0 

NOZZLE  DEFLECTION  (IF  VARIABLE)  7  O 

ATMOSPHERE  .  DEVIATION  iCAO.o 

UUTPUT  l  1-TEXT.  7-NQi..  .  3-fcOTH)  7  3 

0/P  FILENAME  FOR  UNIT  £,  EXMMPL.TXT 

u/p  FILENAME  for  UNIT  Ji  7  EXAMPL.NUM 

calculation  » 

ALTITUDE  HF*  -  0.0  FT 

ALTITUDE  HP  -  1UUUO.O  FT 

STOP 

END  OF  EXECUTION 

CPU  TIME*  7.02  ELAPSED  TIME*  2*34.24 
LXI7 

. RUN  P2V2 

COME’. AT  PERFORMANCE  PROCESSING 
OPTION  OR  (CR)  PGR  HELP  •  42 

42  OMEGA/ETA.  PS  VS  MACH  FOR  GIVEN  HEIGHT.  (P2.PLT) 
DATA  BASE  FILENAME  *  EXAMPL.NUM 

«»*  DATA  ARE  IN  IMPERIAL  UNITS.  ENERGY  PARAMETER  IS  PS 
*««  CONSTANT  INC.  OF  LOAD  FACTOR  *** 

REPLY  •YES’*.  "NO”.  "ALL"  OR  -END"  ' 

HEIGHT  -  O.U  FT  7  N 

HEIGHT  -  10000.0  FT  7  ALL 

END  OF  EXECUTIUN 

CPU  TIME*  5.  OS  ELAPSED  TIME*  5D.04 
EXIT 


FIG.  10.21  Sample  Dialogue  for  Option  4Z 
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APPENDIX  A 


SUMMARY  OF  MODIFICATIONS 

The  modifications  incorporated  are  labelled  Kl,  K2 . K8  in  the 

program  listing  and  are  as  follows: 

1.  MODIFICATION  Kl 

Addition  of  an  option  (AG)  in  program  P2V2  to  plot  contours  of  the  energy 
parameter  ;»  on  a  Mach  number,  turn  rate  (or  load  factor)  grid;  one  page  of  plots 
being  produced  for  each  pressure  height  (or  energy  state). 

This  involved  some  modifications  to  the  routines  in  P2V2  and  two  new 
subroutines  GRIDC  and  PuPSAG  in  PAV2. 

2.  MODIFICATION  K2 

Calculation  and  plotting  of  the  boundaries  corresponding  to  the  lift 
limit  and  structural  limit  in  the  options  AA  and  AG  available  in  P2V2. 

This  involved  some  modifications  in  P1V2  to  calculate  and  print  the 
boundary  points,  and  in  P2V2  to  calculate  the  points  on  the  boundarv  plot  for 
option  AA.  In  addition  two  new  subroutines  BOUND  and  ESTBOUND  where  provided  in 
PAV2  . 

3.  MODIFICATION  K3 

Restructuring  of  subprogram  P1V2  and  the  user  supplied  routines  to 
facilitate  the  programming  of  these  user  routines. 

Also  refers  to  modifications  in  P1V2  and  P2ALIB  to  allow  for  vectored 
thrust  aircraft. 

A.  MODIFICATION  KA 

Provision  of  a  turn  rate  grid  as  an  additional  option  (as  an  alternative 
to  the  load  factor  grid)  when  running  ACRAFT  and  hence  the  option  of  plotting 
either  load  factor  or  turn  rate  when  using  various  plotting  options  available  in 
P2V2. 

5.  MODIFICATION  K3 

The  additional  option  of  a  fourth  type  of  energy  parameter  , „  =  Fc. 'V 
has  involved  modifications  in  P1V2  and  in  P2ALIB. 

6.  MODIFICATION  K6 

Modifications  to  allow  for  a  PACAM2  tabulated  data  file  in  ASCII  to  be 
used  as  input  by  the  user  defined  routines  THRUST  and  AERO  as  an  alternative  to 
the  spline  fitted  information  stored  in  binary. 


i 

1 


I 


7.  MODIFICATION  K7 


Restructuring  of  P4V2  which  was  mostly  rewritten.  Subroutines  PCON,  PES, 
GRIDC,  P4PS4A,  ?4PSiG  are  new  and  other  modifications  were  made  in  P4V2  and 
P4LIB. 

8.  MODIFICATION  K8 

Following  on  from  modification  K4,  P1V2  was  modified  and  RATE2  (in  P 2 V 2 ) 
was  rewritten  in  order  to  improve  the  accuracy  of  the  results  in  option  4C. 
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